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ABSTRACT 

In the first part of this thesis the relativistic viscous fluid equations describing the 
outflow of high temperature matter created via Hawking radiation from microscopic 
black holes are solved numerically for a realistic equation of state. We focus on black 
holes with initial temperatures greater than 100 GeV and lifetimes less than 6 days. 
The spectra of direct photons and photons from ir° decay are calculated for energies 
greater than 1 GeV. We calculate the diffuse gamma ray spectrum from black holes 
distributed in our galactic halo. However, the most promising route for their observa- 
tion is to search for point sources emitting gamma rays of ever-increasing energy. We 
also calculate the spectra of all three flavors of neutrinos arising from direct emission 
from the fluid at the neutrino-sphere and from the decay of pions and muons from 
their decoupling at much larger radii and smaller temperatures for neutrino energies 
between 1 GeV and the Planck energy. The results for neutrino spectra may be ap- 
plicable for the last few hours and minutes of the lifetime of a microscopic black hole. 
In the second part of this thesis the combined field equations of gravity and a scalar 
field are studied. When a potential for a scalar field has two local minima there arise 
spherical shell-type solutions of the classical field equations due to gravitational at- 
traction. We establish such solutions numerically in a space which is asymptotically 
de Sitter. It generically arises when the energy scale characterizing the scalar field 
potential is much less than the Planck scale. It is shown that the mirror image of the 
shell appears in the other half of the Penrose diagram. The configuration is smooth 
everywhere with no physical singularity. 

i 



Acknowledgements 



I would like to express my deep appreciation to Joseph Kapusta for his guidance 
and insight throughout this research. With his kindness and patience, he surpassed 
the role of adviser. He is and will be a role model and a mentor throughout my 
life. I would also like to express my deep appreciation to Yutaka Hosotani for his 
contributions in this thesis. Finally, I would like to thank my wife, Paula, for her 
patience during the years of graduate school. 



ii 



Dedication 

To Paula, Nader, Eghbaleh and Hamoon. 



iii 



Table of Contents 



Abstract i 

Acknowledgements ii 

Dedication iii 

List of Figures vi 

1 Introduction 1 

1.1 Microscopic Black Holes 3 

1.2 Cosmic Shells 6 

2 High Temperature Matter and Particle Spectra from Microscopic 
Black Holes 9 

2.1 Hawking Radiation 9 

2.2 Relativistic Viscous Fluid Equations 11 

2.3 Equation of State and Transport Coefficients 13 

2.4 Numerical Solution and Scaling 16 

2.5 Onset of Free-Streaming 23 

2.6 Photon Emission 28 

2.6.1 Direct photons 28 

2.6.2 7T° decay photons 29 

2.6.3 Instantaneous and integrated photon spectra 30 



iv 



2.7 Observability of Gamma Rays 32 

2.7.1 Diffuse spectra from the galactic halo 32 

2.7.2 Point source systematics 35 

2.8 Neutrino Emission 38 

2.8.1 Directly Emitted Neutrinos 39 

2.8.2 Neutrinos from Pion Decay 43 

2.8.3 Neutrinos from Muon Decay 44 

2.9 Comparison of Neutrino Sources 47 

2.10 Observability of the Neutrino Flux 53 

3 Cosmic Shells 56 

3.1 Shell Solutions in Static Coordinates 56 

3.2 Numerical Analysis of the Nonlinear Regime 66 

3.3 Global Structure of the Space-Time and Solutions 70 

3.4 False Vacuum Decay 78 

4 Summary 80 
Appendix A. Numerical Calculations 84 
References 88 



v 



List of Figures 



2.1 Entropy density as a function of temperature, excluding neutrinos and 
gravitons. It is assumed that the QCD phase transition is first order 

and the EW phase transition is second order 14 

2.2 The value of u — wy at one-tenth the Schwarzschild radius as deter- 
mined by numerical solution. The physical applicability of the numer- 
ical solution begins at radii greater than r s 18 

2.3 The radial dependence of u for three different Hawking temperatures. 
The curves begin at r^/lO and terminate when the local temperature 
reaches 10 MeV 19 

2.4 The radial dependence of T for three different Hawking temperatures. 
The curves begin at r^/lO and terminate when the local temperature 
reaches 10 MeV 20 

2.5 The radial dependence of s for three different Hawking temperatures. 
The curves begin at rs/10 and terminate when the local temperature 
reaches 10 MeV 21 

2.6 The radial dependence of the entropy flow for three different Hawking 
temperatures. The curves begin at r^/lO and terminate when the local 
temperature reaches 10 MeV 22 

2.7 Free-streaming of particles at freeze-out surface 23 

2.8 The rate for mr scattering and for the bremsstrahlung reaction ee — > 
ee7 are compared to the local volume expansion rate. The Hawking 
temperature is 10 TeV 24 



vi 



2.9 The freeze-out or free- streaming radius as a function of the Schwarzschild 
radius for three different freeze-out temperatures 27 

2.10 The instantaneous gamma ray spectrum for three different Hawking 
temperatures assuming Tf = 140 MeV 31 

2.11 Same as figure 9 but with T f = 100 MeV 32 

2.12 The time- integrated gamma ray spectrum starting from the indicated 
Hawking temperature. Here Tf = 140 MeV 33 

2.13 The flux of diffuse gamma rays coming from our galactic halo. The 
normalization is po = 1 pc~ 3 yr _1 . The halo is assumed to be spheri- 
cally symmetric, q — 1; the results for a flattened halo with q = 2 are 
very similar 35 

2.14 The average gamma ray energy as a function of the remaining lifetime 
of the black hole. The times spanned correspond to approximately 400 
GeV <T H < 200 TeV 37 

2.15 The instantaneous direct neutrino spectra emerging from the fluid with 
neutrino-sphere located at T v = 100 GeV compared to the direct Hawk- 
ing radiation. Also shown are neutrinos arising from pion, direct muon 
and indirect muon decays at a decoupling temperature of Tf = 100 
MeV. Here the black hole temperature is T H — 1 TeV. All curves are 

for one flavor of neutrino or anti-neutrino 48 

2.16 Same as figure 2.15 but with T H = 1 TeV 49 

2.17 The instantaneous direct neutrino spectra emerging from the fluid with 
neutrino-sphere located at T v = 100 GeV compared to the direct Hawk- 
ing radiation. Both curves are for one flavor of neutrino or anti-neutrino. 50 

2.18 Same as figure 2.17 but with T H — 110 GeV 51 

2.19 Same as figure 2.17 but with T H = 120 GeV 52 



vii 



2.20 The time integrated neutrino spectra emerging from a microscopic 
black hole. Here the calculation begins when the black hole tempera- 
ture is T = 1 TeV. Either the direct Hawking radiation of neutrinos or 
the direct neutrino emission from a neutrino-sphere at a temperature 
of 100 GeV should be used. All curves are for one species of neutrino 

or anti-neutrino 53 

2.21 Same as figure 2.20 but with T = 10 TeV 54 

3.1 The potential V[<j>] 57 

3.2 Schematic behavior of 0(r) 61 

3.3 Schematic behavior of H(r) 62 

3.4 0(r) of a solution with f/m P = 1.40 x 10~ 3 , Af/f = 0.002, A = 0.01. 
and r are in the units of m P and l P , respectively. The maximum 
value of <p is smaller than f 2 ~ 0.0014 • m P 64 

3.5 H{r) of a solution with f/m P = 1.40 x lO" 3 , Af/f = 0.002, A = 0.01. 
H(r) decreases in two steps in the shell region. H(r) in the de Sitter 
space (= 1 — r 2 /a 2 ) is plotted in a dotted line 66 

3.6 p(r) of a solution with f/m P = 1.40 x 10~ 3 , Af/f = 0.002, A = 0.01. 

p(r) decreases in two steps, from 1 to 0.6256, in the shell region. ... 67 

3.7 The energy density T 00 for a shell solution with f/m P = 1.40 x 10~ 3 , 
Af/f = 0.002, A = 0.01. Too and r are in the units of mp and l P , 
respectively. The energy density is localized in the two shells over the 

de Sitter background 68 

3.8 The radial pressure T n for a shell solution with f/m P = 1.40 x 10~ 3 , 
Af/f = 0.002, A = 0.01. Tn and r are in the units of mp and l P , 
respectively. The pressure is higher and constant between the two shells. 69 

3.9 A solution with double bounces with f/m P = 1.40 x 10~ 3 ,A/// = 
0.002, A = 0.01. A solution with multiple bounces appears at a smaller 
shell radius than a solution with a single bounce 70 



viii 



3.10 The /-dependence of the shell solutions. A = 0.01 and Af/f = 0.002 
are fixed. R, r#, and a are the radius of the shell, the horizon length 
where H(r) vanishes, and the horizon length in the de Sitter space, 
respectively. The ratios of various pairs are plotted. Circles, squares, 

and diamonds correspond to data points obtained 71 

3.11 Penrose diagram of the de Sitter space 74 

3.12 Global structure of the shell solution in the conformal coordinates. The 
mirror image appears in the other hemisphere 77 

A.l 47rr 2 T 0r /L as a function of r/r s for a black hole temperature of 10 
TeV. The curve begins at r s and terminates when the local temperature 
reaches T = m e 87 



ix 



Chapter 1 
Introduction 



Black holes can be created in the collapse of a star which was originally 10 to 
20 times the mass of the sun. There is strong evidence that there are huge black 
holes of about 1 million solar masses near the center of galaxies. Another possibility 
are primordial black holes of almost arbitrarily small mass which were created in the 
early universe. Depending on the time of the formation of primordial black holes in 
the early universe, they can have a wide range of masses. The black holes which 
were formed at the the Planck time (10~ 43 s) would have the Planck mass (10~ 5 
g), and those were formed at 1 s would be as massive as 10 5 solar masses. It was 
shown by Hawking in 1974 that a black hole radiates thermally with a temperature 
that is inversely proportional to its mass. This phenomenon, which is called Hawk- 
ing radiation, becomes important for black holes with small enough mass. In fact, 
primordial black holes with a mass less than about 10 15 g have already evaporated 
within the lifetime of the universe, which is 10 10 yrs. The size of such black holes is 
less than 10~ 16 m, making them microscopic. Primordial black holes with an initial 
mass a little more than 10 15 g are evaporating at the present time at a rate which 
makes it possible to detect them. In the first part of this thesis we study the particle 
spectra emitted from primordial microscopic black holes in the last 6 days of their 
existence. We assume that at very high black hole temperatures the particles emitted 
from black hole via Hawking radiation will scatter from each other enough to create 
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a shell of thermalized matter surrounding these black holes. This shell of hot mat- 
ter will change the emergent spectrum from microscopic black holes. More particles 
with lower average energy will be emitted out of the thermalized matter compared 
to fewer particles with higher average energy emitted directly via Hawking radiation. 
This phenomenon is similar to the photosphere of stars. Because of the photosphere 
the temperature on the outer surface of a star is much lower than the temperature 
at the core of the star. In this thesis we show that the most promising route for the 
observation of microscopic black holes is to search for point sources emitting gamma 
rays of ever-increasing energy. A black hole with a temperature above 100 GeV and a 
Schwarzschild radius less than 10~ 4 fm will get brighter on a time scale of 6 days and 
then disappear. Such an observation would be remarkable, possibly unique, because 
astrophysical sources of gamma rays normally cool at late times. This would directly 
reflect the increasing Hawking temperature as the black hole explodes and disap- 
pears. Since the highest temperatures in the present universe exist in the vicinity of 
microscopic black holes, by studying the emergent spectrum from these black holes 
we can learn about the physics at very high temperatures which is not achievable 
in terrestrial experiments. Observation of microscopic black holes will shed light on 
the biggest problems in physics today. By studying and observing microscopic black 
holes we will be able to test the physics above the electroweak scale and beyond the 
four dimensional Standard Model. 

It is possible that the source of the dark energy in our universe is a scalar field. 
Therefore, it is important to study the solutions to the combined field equations of 
gravity and a scalar field. In the second part of this thesis we report the discovery of 
static shell-like solutions to the combined field equations of gravity and a scalar field 
with a potential which has two non-degenerate minima. The absolute minimum of the 
potential is the true vacuum and the other minimum is the false vacuum. Both inside 
and outside of the shells are de Sitter space and their radii are comparable to the size 
of the horizon of the universe. The shells are spherical and very thin compared to 
their radii. The energy density in the shells is much higher than the energy density of 



3 

the de Sitter space inside and outside of the shells. If anything like these structures 
exist in nature they most likely would have been created in the early universe and 
are therefore cosmological. We know of no other way to produce them. 

In chapter 2 of this thesis we present the work done in collaboration with J. 
Kapusta which can be found in [1] and [2]. In chapter 3 we present the work done in 
collaboration with Y. Hosotani, J. Kapusta and T. Nakajima which can be found in 
[3]. 

1.1 Microscopic Black Holes 

Hawking radiation from black holes [4] is of fundamental interest because it is a 
phenomenon arising out of the union of general relativity and quantum mechanics, 
a situation that could potentially be observed. It is also of great interest because 
of the temperatures involved. A black hole with mass M radiates thermally with a 
Hawking temperature T H = ml/8irM where m P = G~ 1/2 = 1.22 x 10 19 GeV is the 
Planck mass. (Units are k — c — k B — 1.) In order for the black hole to evaporate 
rather than accrete it must have a temperature greater than that of the present-day 
blackbody radiation of the universe of 2.7 K = 2.3 xlO~ 4 eV. This implies that M 
must be less than 1% of the mass of the Earth. 

Such small black holes most likely would have been formed primordially; there is no 
other mechanism known to form them. As the black hole radiates, its mass decreases 
and its temperature increases until T# becomes comparable to the Planck mass, at 
which point the semiclassical calculation breaks down and the regime of full quantum 
gravity is entered. Only in two other situations are such enormous temperatures 
achievable: in the early universe and in central collisions of heavy nuclei like gold or 
lead. Even then only about T = 500 MeV is reached at the RHIC (Relativistic Heavy 
Ion Collider) just completed at Brookhaven National Laboratory and T = 1 GeV is 
expected at the LHC (Large Hadron Collider) at CERN to be completed in 2006. 
Supernovae and newly formed neutron stars only reach temperatures of a few tens 
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of MeV. To set the scale from fundamental physics, we note that the spontaneously 
broken chiral symmetry of QCD gets restored in a phase transition or rapid crossover 
at a temperature around 170 MeV, while the spontaneously broken gauge symmetry 
in the electroweak sector of the standard model gets restored in a phase transition or 
rapid crossover at a temperature around 100 GeV. The fact that temperatures of the 
latter order of magnitude will never be achieved in a terrestrial experiment should 
motivate us to study the fate of microscopic black holes during the final days, hours 
and minutes of their lives when their temperatures have risen to 100 GeV and above. 
In this thesis we shall focus on Hawking temperatures greater than 100 GeV. The 
fact that microscopic black holes have not yet been observed [5] should not be viewed 
as a deterrent, but rather as a challenge for the new millennium! 

There is some uncertainty over whether the particles scatter from each other after 
being emitted, perhaps even enough to allow a fluid description of the wind coming 
from the black hole. Let us examine what might happen as the black hole mass 
decreases and the associated Hawking temperature increases. 

When Th <C rn e (electron mass) only photons, gravitons, and neutrinos will be 
created with any significant probability. These particles will not significantly interact 
with each other, and will freely propagate away from the black hole with energies of 
order T H . Even when T H m m c the Thomson cross section is too small to allow the 
photons to scatter very frequently in the rarified electron-positron plasma around the 
black hole. This may change when Th ~ 100 MeV when muons and charged pions 
are created in abundance. At somewhat higher temperatures hadrons are copiously 
produced and local thermal equilibrium may be achieved, although exactly how is 
an unsettled issue. Are hadrons emitted directly by the black hole? If so, they will 
be quite abundant at temperatures of order 150 MeV because their mass spectrum 
rises exponentially (Hagedorn growth as seen in the Particle Data Group tables [6]). 
Because they are so massive they move nonrelativistically and may form a very dense 
equilibrated gas around the black hole. But hadrons are composites of quarks and 
gluons, so perhaps quark and gluon jets are emitted instead? These jets must decay 
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into the observable hadrons on a typical proper length scale of 1 fm and a typical 
proper time scale of 1 fm/c. This was first studied by MacGibbon and Webber 
[7] and MacGibbon and Carr [8]. Subsequently Heckler [9] argued that since the 
emitted quarks and gluons are so densely packed outside the event horizon they are 
not actually fragmenting into hadrons in vacuum but in something more like a quark- 
gluon plasma, so perhaps they thermalize. He also argued that QED bremsstrahlung 
and pair production were sufficient to lead to a thermalized QED plasma when Th 
exceeded 45 GeV [10]. These results are somewhat controversial and need to be 
confirmed. The issue really is how to describe the emission of wave packets via the 
Hawking mechanism when the emitted particles are (potentially) close enough to be 
mutually interacting. A more quantitative treatment of the particle interactions on 
a semiclassical level was carried out by Cline, Mostoslavsky and Servant [11]. They 
solved the relativistic Boltzmann equation with QCD and QED interactions in the 
relaxation-time approximation. It was found that significant particle scattering would 
lead to a photosphere though not perfect fluid flow. 

Rather than pursuing the Boltzmann transport equation Kapusta applied rela- 
tivistic viscous fluid equations to the problem assuming sufficient particle interaction 
[12]. It was found that a self-consistent description emerges of a fluid just marginally 
kept in local thermal equilibrium, and that viscosity is a crucial element of the dynam- 
ics. The fluid description has been used to calculate the particle spectra emitted in 
heavy ion collisions and has been very successful in this application [13]. The purpose 
of this thesis is a more extensive analysis of these equations and their observational 
consequences. 

The plan of the second chapter of this thesis is as follows. In Sec. 2.1 we give a 
brief review of Hawking radiation sufficient for our uses. In Sec. 2.2 we give the set of 
relativistic viscous fluid equations necessary for this problem along with the assump- 
tions that go into them. In Sec. 2.3 we suggest a relatively simple parametrization 
of the equation of state for temperatures ranging from several MeV to well over 100 
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GeV. We also suggest a corresponding parametrization of the bulk and shear vis- 
cosites. In Sec. 2.4 we solve the equations numerically, study the scaling behavior 
of the solutions, and check their physical self-consistency. In Sec. 2.5 we estimate 
where the transition from viscous fluid flow to free-streaming takes place. In Sec. 
2.6 we calculate the instantaneous and time-integrated spectra of high energy pho- 
tons from the two dominant sources: direct and neutral pion decay. In Sec. 2.7 we 
study the diffuse gamma ray spectrum from microscopic black holes distributed in 
our galactic halo. We also study the systematics of gamma rays from an individual 
black hole, should we be so fortunate to observe one. In Sec. 2.8 we calculate the 
instantaneous and time-integrated spectra of high energy neutrinos from the three 
dominant sources: direct, pion decay and muon decay. In Sec. 2.9 we compare the 
spectra from all neutrino sources graphically. We also compare our results with the 
spectra of neutrinos emitted directly as Hawking radiation without any subsequent 
interactions. In Sec. 2.10 we study the possibility of observing neutrinos from a 
microscopic black hole directly. 

1.2 Cosmic Shells 

Gravitational interactions, which are inherently attractive for ordinary matter, can 
produce soliton-like objects even when they are strictly forbidden in flat space. They 
are possible as a consequence of the balance between repulsive and attractive forces. 
One such example is a monopole or dyon solution in the pure Einstein- Yang-Mills 
theory in the asymptotically anti-de Sitter space [14, 15, 16]. In the pure Yang-Mills 
theory in flat space there can be no static solution at all [17] but once gravitational 
interaction is included there arise particle like solutions [18]. Whereas all solutions 
are unstable in the asymptotically flat or de Sitter space, there appear a continuum 
of stable monopole and dyon solutions in the asymptotically anti-de Sitter space. 
The stable solutions are cosmological in nature; their size is typically of order |A| -1 / 2 
where A is the cosmological constant. 
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The possibility of false vacuum black holes has also been explored. Suppose that 
the potential in a scalar field theory has two minima, one corresponding to the true 
vacuum and the other to the false vacuum. If the universe is in the false vacuum, 
a bubble of the true vacuum is created by quantum tunneling which expands with 
accelerated velocity. The configuration is called a bounce [19]. Now flip the configu- 
ration [20]. The universe is in the true vacuum with potential V = and the inside 
of a sphere is excited to the false vacuum with V > 0. Is such a de Sitter lump 
in Minkowski space possible? If the lump is too small it would be totally unstable. 
The energy localized inside the lump can dissipate to spatial infinity. If the lump 
is big enough the Schwarzschild radius becomes larger than the lump radius so that 
the lump is inside a black hole. The energy cannot escape to infinity. It looks like 
a soliton in Minkowski space. However, as a black hole it is a dynamical object. 
The configuration is essentially time dependent. This false vacuum black hole con- 
figuration, however, does not solve the static field equations at the horizon. It has 
been recently proven that there can be no such everywhere-regular black hole solution 
[21, 22]. Rather, false vacuum lumps in flat space evolve dynamically [23, 24, 25]. 

The purpose of this part of thesis is to report new solutions to the coupled equa- 
tions of gravity and scalar field theory which display a spherically symmetric shell 
structure [26]. We demonstrate that such a structure appears when the potential for 
a scalar field has two local minima and the space is asymptotically de Sitter. In the 
examples we present, both the inside and outside of the shells are de Sitter space 
with the same cosmo logical constant. The structure becomes possible only when the 
energy scale of the scalar field potential becomes small compared with the Planck 
scale. While such shell structures might not be easy to create in the present universe, 
it is quite plausible that they could have been created during a phase transition early 
in the universe [23]. A similar configuration has been investigated in Ref. [27]. 

The plan of the third chapter of this thesis is as follows. In Sec. 3.1 we precisely 
state the problem, solve the field equations in those regions of space-time where 
they can be linearized, and sketch the solution in the nonlinear shell region in static 
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coordinates. In Sec. 3.2 we solve the nonlinear equations in the shell region and 
display the dependence on the parameters of the theory. In Sec. 3.3 we extend 
the solution from static coordinates, which have a coordinate singularity, to other 
coordinate systems that do not, thereby displaying the existence and character of the 
solution throughout the full space-time manifold. In Sec. 3.4 we study the stability of 
the classical solution to quantum fluctuations. Finally in chapter 4 we will conclude 
and summarize both parts of this thesis. 



Chapter 2 



High Temperature Matter and 
Particle Spectra from Microscopic 
Black Holes 



2.1 Hawking Radiation 

There are at least two intuitive ways to think about Hawking radiation from black 
holes. One way is vacuum polarization. Particle-antiparticle pairs are continually 
popping in and out of the vacuum, usually with no observable effect. In the presence 
of matter, however, their effects can be observed. This is the origin of the Lamb effect 
first measured in atomic hydrogen in 1947. When pairs pop out of the vacuum near 
the event horizon of a black hole one of them may be captured by the black hole and 
the other by necessity of conservation laws will escape to infinity with positive energy. 
The black hole therefore has lost energy - it radiates. Due to the general principles 
of thermodynamics applied to black holes it is quite natural that it should radiate 
thermally. An intuitive argument that is more quantitative is based on the uncertainty 
principle. Suppose that we wish to confine a massless particle to the vicinity of a black 
hole. Given that the average momentum of a massless particle at temperature T is 



9 



10 

approximately nT, the uncertainty principle requires that confinement to a region the 
size of the Schwarzschild diameter places a restriction on the minimum value of the 
temperature 

7rT-2r 5 >l/2. (2.1) 

The minimum is actually attained for the Hawking temperature. The various physical 
quantities are related as r s = 2M/mp = 1/A.ttTh- 

The number of particles of spin s emitted with energy E per unit time is given by 
the formula 

dN s T s 1 
dEdt ~ 27 exp(E/T H ) - (-l) 2s ' 

All the computational effort really goes into calculating the absorption coefficient T s 
from a relativistic wave equation in the presence of a black hole. Integrating over all 
particle species yields the luminosity 

L = ~ = a(M)^ 2 = ^a(Tu)Tl . (2.3) 

Here ct(M) is a function reflecting the species of particles available for creation in the 
gravitational field of the black hole. It is generally sufficient to consider only those 
particles with mass less than T H ; more massive particles are exponentially suppressed 
by the Boltzmann factor. Then 



a = 2.011 x 10~ 8 



3700iVo + 2035Ari/ 2 + 835Aq + 95iV 2 (2.4) 



where N s is the net number of polarization degrees of freedom for all particles with 
spin s and with mass less than T#. The coefficients for spin 1/2, 1 and 2 were 
computed by Page [28] and for spin by Sanchez [29]*. In the standard model 
N — 4 (Higgs boson), N 1 / 2 = 90 (three generations of quarks and leptons), Aq = 
24 [SU(3)xSU(2)xU(l) gauge theory], and N 2 = 2 (gravitons). This assumes T H 
is greater than the temperature for the electroweak gauge symmetry restoration. 



*In [1] the coefficient for spin particles was taken as 4200. In Eq. (2.4) we have reduced this 
value to 3700 which results in a small correction of 0.7% in the numerical value of <x(Th > 100 GeV). 
We thank Harald Anlauf for providing us with this more accurate coefficient 
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Numerically a (T H > 100 GeV) =4.4xl0~ 3 . Starting with a black hole of temperature 
T H , the time it takes to evaporate or explode is 

2 

Trip 

At = 3a(T H )(87rT H y ' (2 ' 5) 

This is also the characteristic time scale for the rate of change of the luminosity of a 
black hole with temperature T#. 

At present a black hole will explode if T# > 2.7 K and correspondingly M < 
4.6 x 10 25 g which is approximately 1% of the mass of the Earth. More massive 
black holes are cooler and therefore will absorb more matter and radiation than they 
radiate, hence grow with time. Taking into account emission of gravitons, photons, 
and neutrinos a critical mass black hole today has a Schwarszchild radius of 68 microns 
and a lifetime of 2 x 10 43 years. 



2.2 Relativist ic Viscous Fluid Equations 

The relativistic imperfect fluid equations describing a steady-state, spherically 
symmetric flow with no net baryon number or electric charge and neglecting gravity 
(see below) are T^ v = black hole source. The nonvanishing components of the energy- 
momentum tensor in radial coordinates are [30] 

T 00 = 1 2 {p + e) _ p + v 2 ATd . ss 

T 0r = V1 2 (P + e)+vAT diss 

T rr = v ^ 2 (p + e) + P + AT diss (2.6) 

representing energy density, radial energy flux, and radial momentum flux, respec- 
tively, in the rest frame of the black hole. Here v is the radial velocity with 7 the 
corresponding Lorentz factor, u = vy, e and P are the local energy density and 
pressure, and 

4 n / du u\ a / du 2u\ . 
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where 77 is the shear viscosity and ( is the bulk viscosity. A thermodynamic identity 
gives Ts = P + e for zero chemical potentials, where T is temperature and s is entropy 
density. There are two independent differential equations of motion to solve for the 
functions T(r) and v(r). These may succinctly written as 

^-(r 2 T 0r ) = 
dr v > 

^-(r 2 T rr ) = 0. (2.8) 
dr v ' 

An integral form of these equations is sometimes more useful since it can readily 
incorporate the input luminosity from the black hole. The first represents the equality 
of the energy flux passing through a sphere of radius r with the luminosity of the black 
hole: 

47rr 2 T 0r = L ( 2 _g) 

The second follows from integrating a linear combination of the differential equations. 
It represents the combined effects of the entropy from the black hole together with 
the increase of entropy due to viscosity: 

rr \ 

Aitr 2 us = Aix I dr'r' 2 — 

Jro T 

The term L/Th arises from equating the entropy per unit time lost by the black 
hole —dSbh/dt with that flowing into the matter. Using the area formula for the 
entropy of a black hole, = rripTcrg = 4itM 2 /mp, and identifying —dM/dt with 
the luminosity, the entropy input from the black hole is obtained. 

The above pair of equations are to be applied beginning at some radius r greater 
than the Schwarzschild radius r s , that is, outside the quantum particle production 
region of the black hole. The radius r at which the imperfect fluid equations are first 
applied should be chosen to be greater than the Schwarzschild radius, otherwise the 
computation of particle creation by the black hole would be invalid. It should not 
be too much greater, otherwise particle collisions would create more entropy than is 
accounted for by the equation above. The energy and entropy flux into the fluid come 



8 (du u \ 2 ^ / du ^ 2u\ 2 

9 ^ V dr' r' 1 \ dr' r' ) 



L 



(2.10) 
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from quantum particle creation by the black hole at temperature T H . Gravitational 
effects are of order r s /r, hence negligible for r > (5 — I0)r s . 



2.3 Equation of State and Transport Coefficients 

Determination of the equation of state as well as the two viscosities for tem- 
peratures ranging from MeV to TeV and more is a formidable task. Here we shall 
present some relatively simple parametrizations that seem to contain the essential 
physics. Improvements to these can certainly be made, but probably will not change 
the viscous fluid flow or the observational consequences very much. 

The hot shell of matter surrounding a primordial black hole provides a theoretical 
testing ground rivaled only by the big bang itself. To illustrate this we have plotted 
a semi-realistic parametrization of the equation of state in figure 2.1. Gravitons and 
neutrinos are not included. We assume a second order electroweak phase transition at 
a temperature of T E w = 100 GeV. Above that temperature the standard model has 
101.5 effective massless bosonic degrees of freedom (as usual fermions count as 7/8 of 
a boson). We assume a first order QCD phase transition at a temperature of Tqcd = 
170 MeV. The number of effective massless bosonic degrees of freedom changes from 
47.5 just above this critical temperature (u, d, s quarks and gluons) to 7.5 just below it 
(representing the effects of all the massive hadrons in the Particle Data Group tables) 
[31]. Below 30 MeV only electrons, positrons, and photons remain, and finally below 
a few hundred keV only photons survive in any appreciable number. The explicit 
parametrization shown in figure 2.1 is as follows: 



4tt 2 

s m = —t 3 i 

V ; 90 



101.5, Tew < T, 

56.5 + 45e-( T — " T )/ T , T QCD <T< T EW , (2.11) 

2 + 3.5 e~ m ^ T + 27.25 <tVqcd-T)/t ? T < Tqcd 

It may very well be that there are no true thermodynamic phase transitions in the 
standard model but only rapid crossovers from one phase to the other. None of our 
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Figure 2.1: Entropy density as a function of temperature, excluding neutrinos and 
gravitons. It is assumed that the QCD phase transition is first order and the EW 
phase transition is second order. 

calculations or results depend on such details. A word about neutrinos: It is quite 
possible that they should be considered in approximate equilibrium at temperatures 
above 100 GeV where the electroweak symmetry is restored. Still there is some 
uncertainty about this. Since they provide only a few effective degrees of freedom out 
of more than 100 their neglect should cause negligible error. 

Now we turn to the viscosities. The shear viscosity was calculated in [32] for the 
full standard model in the symmetry restored phase, meaning temperatures above 
100 GeV or so, using the relaxation time approximation. The result is 

r](T > 100 GeV) = 82.5 T 3 (2.12) 



15 



when numerical values for coupling constants etc. are put in. The shear viscosity for 
QCD degrees of freedom only was calculated to leading order in the QCD coupling 
a s in [33] to be 



where Nf is the number of quark flavors whose mass is less than T. An improved 
calculation for gauge theories was given in [34] ; for QCD there is very little difference 
with [33]. We observe that the ratio of the shear viscosity to the entropy density 
as appropriate for the above two cases, is dimensionless and has about the same 
numerical value in both. Therefore, as a practical matter we assume that the shear 
viscosity always scales with the entropy density for all temperatures of interest. We 
take the constant of proportionality from the full standard model cited above: 



There is even less known about the bulk viscosity at the temperatures of interest to 
us. The bulk viscosity is zero for point particles with no internal degrees of freedom 
and with local interactions among them. In renormalizable quantum field theories 
the interactions are not strictly local. In particular, the coupling constants acquire 
temperature dependence according to the renormalization group. For example, to one 
loop order the QCD coupling has the functional dependence a s ~ l/ln(T/A) where 
A is the QCD scale. On account of this dependence the bulk viscosity is nonzero. We 
estimate that 



and this is what we shall use in the numerics. 

Overall we have a modestly realistic description of the equation of state and the 
viscosities that are still a matter of theoretical uncertainty. One needs s(T) : rj(T) : ((T) 
over a huge range of T. Of course, these are some of the quantities one hopes to obtain 
experimental information on from observations of exploding black holes. 




(2.13) 




(2.14) 



C « 1(T 4 77 



(2.15) 
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2.4 Numerical Solution and Scaling 

Several limiting cases of the relativistic viscous fluid equations were studied in [12]. 
The most realistic situation used the equation of state e = aT 4 , s = (4/3)aT 3 and 
viscosities 77 = 65T 3 , ( = &bT 3 with the coefficients a, bs, &b all constant. A scaling 
solution, valid at large radii when 7 ^> 1, was found to be T(r) = T (r /r) 2 / 3 and 
7(7-) = 7o(?"/Vo) 1 ^ 3 - The constants must be related by 36aT r = (326 s + 441&b)7o- 
This r-dependence of T and 7 is exactly what was conjectured by Heckler [10]. 

It was shown in [12] that if approximate local thermal equilibrium is achieved 
it can be maintained, at least for the semi-realistic situation described above. The 
requirement is that the inverse of the local volume expansion rate 9 = -u M .^ be compa- 
rable to or greater than the relaxation time for thermal equilibrium [30]. Expressed 
in terms of a local volume element V and proper time r it is 9 = (l/V)dV/dr, 
whereas in the rest frame of the black hole the same quantity can be expressed as 
(l/r 2 )d(r 2 u)/dr. Explicitly 



Of prime importance in achieving and maintaining local thermal equilibrium in a 
relativistic plasma are multi-body processes such as 2 — > 3 and 3^2, etc. This 
has been well-known when calculating quark-gluon plasma formation and evolution 
in high energy heavy ion collisions [35, 36] and has been emphasized in Refs. [9, 10] in 
the context of black hole evaporation. This is a formidable task in the standard model 
with its 16 species of particles. Instead three estimates for the requirement that local 
thermal equilibrium be maintained were made. The first and simplest estimate is to 
require that the thermal de Broglie wavelength of a massless particle, 1/3T, be less 
than 1/9. The second estimate is to require that the Debye screening length for each of 
the gauge groups in the standard model be less than 1/9. The Debye screening length 
is the inverse of the Debye screening mass where n = 1, 2, 3 for the gauge groups 
U(l), SU(2), SU(3). Generically oc g n T where g n is the gauge coupling constant 
and the coefficient of proportionality is essentially the square root of the number of 




(2.16) 
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charge carriers [37]. For example, for color SU(3) = ^yl + Nf/6T where Nf is 
the number of light quark flavors at the temperature T. The numerical values of the 
gauge couplings are: g\ = 0.344, g 2 = 0.637, and g 3 = 1.18 (evaluated at the scale 
mz) [6]. So within a factor of about 2 we have m D « T. The third and most relevant 
estimate is the mean time between two-body collisions in the standard model for 
temperatures greater than the electroweak symmetry restoration temperature. This 
mean time was calculated in [32] in the process of calculating the viscosity in the 
relaxation time approximation. Averaged over all particle species in the standard 
model one may infer from that paper an average time of 3.7/T. Taking into account 
multi-body reactions would decrease that by about a factor of two to four. All three of 
these estimates are consistent within a factor of 2 or 3. The conclusion to be drawn 
is that local thermal equilibrium should be achieved when 9 < T. Once thermal 
equilibrium is achieved it is not lost because 6/T is independent of r. The picture 
that emerges is that of an imperfect fluid just marginally kept in local equilibrium by 
viscous forces. 

The results quoted above are only valid at large r and for the equation of state 
s cc T 3 . To know the behavior of the solution at non- asymptotic r and for the 
more sophisticated equation of state and viscosities described in Sec. 2.3 requires a 
numerical analysis. We have found that the most convenient form of the viscous fluid 
equations for numerical evaluation are 



for entropy flow [from Eq. (2.10)]. Obviously the entropy flux is a monotonically 
increasing function of r because of dissipation. 

Mathematically the above pair of equations apply for all r > 0, although physically 
we should only apply them beyond the Schwarzschild radius rs- Let us study them 




(2.17) 



for energy conservation [from Eq. (2.9)] and 




(2.18) 



18 

first in the limit r — > 0, which really means the assumption that t)<l. Then u ~ t> 
and 7«1. We also consider black hole temperatures greater than Tew so that the 
equation of state and the viscosities no longer change their functional forms. It is 
straightforward to check that a power solution satisfies the equations, with 



u{r) = Uiir/ri) 2 ^ 
T(r) = T^/rf 5 



(2.19) 



where r\ is some reference radius. If the luminosity and the reference radius are given 
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Figure 2.2: The value of u = v 7 at one-tenth the Schwarzschild radius as determined 
by numerical solution. The physical applicability of the numerical solution begins at 
radii greater than r s . 



then Ui and Tj are determined by the fluid equations. 

The numerical solution for all radii needs some initial conditions. Typically we 
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begin the solution at one-tenth the Schwarzschild radius. At this radius the itj, as 
determined above, is small enough to serve as a good first estimate. However, it needs 
to be fine-tuned to give an acceptable solution at large r. For example, at large r 
there is an approximate but false solution: T = const with u ~ r. The problem is 
that we need a solution valid over many orders of magnitude of r. If Eq. (2.17) is 
divided by r 2 and if the right hand side is neglected in the limit r — > oo then the 
left hand side is forced to be zero. We have used both Mathematica and a fourth 
order Runge-Kutta method with adaptive step-size to solve the equations. They give 
consistent results. For more details on numerical calculations, see appendix A. 




Figure 2.3: The radial dependence of u for three different Hawking temperatures. The 
curves begin at r^/lO and terminate when the local temperature reaches 10 MeV. 
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Figure 2.4: The radial dependence of T for three different Hawking temperatures. 
The curves begin at rg/10 and terminate when the local temperature reaches 10 
MeV. 

In figure 2.2 we plot w(rs/10) versus T# . It is essentially constant for all T# > Tew 
with the value of 0.0415. In figure 2.3 we plot the function u(r) versus r for three 
different black hole temperatures. The radial variable r is expressed in units of its 
value when the temperature first reaches Tqcd, and u is expressed in units of its value 
at that same radius. This allows us to compare different black hole temperatures. To 
rather good accuracy these curves seem to be universal as they essentially lie on top 
of one another. The curves are terminated when the temperature reaches 10 MeV. 
The function u(r) behaves like r 1//3 until temperatures of order 100 MeV are reached. 
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Figure 2.5: The radial dependence of s for three different Hawking temperatures. The 
curves begin at rg/10 and terminate when the local temperature reaches 10 MeV. 

The simple parametrization 

u(r) =u s (r/r s ) 1/3 (2.20) 

with us = 0.10 will be very useful when studying radiation from the surface of the 
fluid. 

In figure 2.4 we plot the temperature in units of Tq C d versus the radius in units 
°f r QCD for the same three black hole temperatures as in figure 2.3. Again the 
curves are terminated when the temperature drops to 10 MeV. The curves almost 
fall on top of one another but not perfectly. The temperature falls slightly slower 
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than the power-law behavior r~ 2 / 3 expected on the basis of the equation of state 
s = (4/3)aT 3 . The reason is that the effective number of degrees of freedom is falling 
with the temperature. The entropy density is shown in figure 2.5. It also exhibits an 
imperfect degree of scaling similar to the temperature. 




Figure 2.6: The radial dependence of the entropy flow for three different Hawking 
temperatures. The curves begin at r^/lO and terminate when the local temperature 
reaches 10 MeV. 

Since viscosity plays such an important role in the outgoing fluid we should expect 
significant entropy production. In figure 2.6 we plot the entropy flow Awr US cLS cL 
function of radius for the same three black hole temperatures as in figures 2.3-2.5. It 
increases by several orders of magnitude. The fluid flow is far from isentropic. 



2.5 Onset of Free- Streaming 



23 



Eventually the fluid expands so rapidly that the particles composing the fluid lose 
thermal contact with each other and begin free-streaming as is shown in figure 2.7. 
In heavy ion physics this is referred to as thermal freeze-out, and in astronomy it is 
usually associated with the photosphere of a star. In the sections above we argued 




Figure 2.7: Free-streaming of particles at freeze-out surface. 



that thermal contact should occur for all particles, with the exception of gravitons 
and neutrinos, down to temperatures on the order of Tqcd- Below that temperature 
the arguments given no longer apply directly; for example, the relevant interactions 



are not those of perturbative QCD. 
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Figure 2.8: The rate for tttt scattering and for the bremsstrahlung reaction ee — > ecy 
are compared to the local volume expansion rate. The Hawking temperature is 10 
TeV. 

Extensive studies have been made of the interactions among hadrons at finite 
temperature. Prakash et al. [38] used experimental information to construct scat- 
tering amplitudes for pions, kaons and nucleons and from them computed thermal 
relaxation rates. The relaxation time for n — n scattering can be read off from their 
figures and simply parametrized as 

T 



w 16 , 

MOO MeV 



MeV 



(2.21) 



which is valid for 100 < T < 200 MeV. This rate is compared to the volume expansion 
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rate 9 (see Sec. 2.4) in figure 2.8. From the figure it is clear that pions cannot maintain 
thermal equilibrium much below 160 MeV or so. Since pions are the lightest hadrons 
and therefore the most abundant at low temperatures, it seems unlikely that other 
hadrons could maintain thermal equilibrium either. 

Heckler has argued vigorously that electrons and photons should continue to in- 
teract down to temperatures on the order of the electron mass [9, 10]. Multi-particle 
reactions are crucial to this analysis. Let us see how it applies to the present situ- 
ation. Consider, for example, the cross section for ee — > eej. The energy-averaged 
cross section is [10] 

<Xbrcm = 8a EM r 2 \n{2E/m e ) (2.22) 

where m e is the electron mass, r = oiEM/m e is the classical electron radius, and 
E ^> m e is the energy of the incoming electrons in the center-of-momentum frame. 
(If one computes the rate for a photon produced with the specific energies 0.1-E 1 , 0.25-E 1 , 
or 0.5-E 1 the cross section would be larger by a factor 4.73, 2.63, or 1.27, respectively.) 
The rate using the energy-averaged cross section is 

3 



-l 

' brem 



2 m T 

7T Z 



8«Im ln(6T/me) 



m 2 



(2.23) 



where we have used the average energy (E) 3T for electrons with m e <C T. This 
rate is also plotted in figure 2.8. It is large enough to maintain local thermal equi- 
librium down to temperatures on the order of 140 MeV. Of course, there are other 
electromagnetic many-particle reactions which would increase the overall rate. On 
the other hand, as pointed out by Heckler [10], these reactions are occurring in a high 
density plasma with the consequence that dispersion relations and interactions are 
renormalized by the medium. If one takes into account only renormalization of the 
electron mass, such that m\^ ~ m\ + e 2 T 2 /3 when m e <C T, then the rate would be 
greatly reduced. 

Does this mean that photons and electrons are not in thermal equilibrium at the 
temperatures we have been discussing? Consider bremsstrahlung reactions in the 
QCD plasma. There are many 2 3 reactions, such as: qiq 2 — > <?i<?27, 
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99 ~ ¥ 997) an d so on - Here the subscripts label the quark flavor, which may or may not 
be the same. The rate for these can be estimated using known QED and QCD cross 
sections [39, 40, 41]. Using an effective quark mass given by gT we find that the rate 
is a s T with a coefficient of order or larger than unity. Since a s becomes of order unity 
near Tq C d we conclude that photons are in equilibrium down to temperatures of that 
order at least. To make the matter even more complicated we must remember that 
the expansion rate 9 is based on a numerical solution of the viscous fluid equations 
which assume a constant proportionality between the shear and bulk viscosities and 
the entropy density. Although these proportionalities may be reasonable in QCD and 
electroweak plasmas at high temperatures they may fail at temperatures below Tqcd- 
The viscosities should be computed using the relaxation times for self-consistency of 
the transition from viscous fluid flow to free-streaming, which we have not done. For 
example, the first estimate for the shear viscosity for massless particles with short 
range interactions is T 4 r where r is the relaxation time. For pions we would get r\ ~ 
const, not r\ ~ T 3 . As another example, we must realize that the bulk viscosity can 
become significant when the particles can be excited internally. This is, in fact, the 
case for hadrons. Pions, kaons and nucleons are all the lowest mass hadrons each 
of which sits at the base of a tower of resonances [6]. See, for example, [42] and 
references therein. 

In order to do gamma ray phenomenology we need a practical criterion for the 
onset of free-streaming. We shall assume that this happens suddenly at a temper- 
ature Tf in the range 100 to 140 MeV. We shall assume that particles whose mass 
is significantly greater than Tf have all annihilated, leaving only photons, electrons, 
muons and pions. In figure 2.9 we plot the freeze-out radius 77 = r(T/) for Tf = 100, 
120 and 140 MeV versus the Schwarzschild radius. The fact that 77 increases as Tf 
decreases is an obvious consequence of energy conservation. More interesting is the 

]/2 1/2 

power-law scaling: 77 ~ r s ~ TjJ . This scaling can be understood as follows. 
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Figure 2.9: The freeze-out or free-streaming radius as a function of the Schwarzschild 
radius for three different freeze-out temperatures. 

The luminosity from the decoupling or freeze-out surface is 

L, = 4vrrJ (^T?7?) d f (2.24) 

where the quantity in parentheses is the surface flux for one massless bosonic degree 
of freedom and df is the total number of effective massless bosonic degrees of freedom. 
For the particles listed above we have df = 12. By energy conservation this is to be 
equated with the Hawking formula for the black hole luminosity, 

L h = Mir 2 a e h ff T 2 H , (2.25) 

where a^f does not include the contribution from gravitons and neutrinos. Together 
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with the scaling function for the flow velocity, Eq. (2.20), we can solve for the radius 

2 /45vraF\ 3/8 [f-, 



H 



and for the boost 



(2.26) 



7/T/ = 2US \2^di J ^ TfTH * °' 22 ^ ^ ' (2 ' 27) 
From these we see that the final radius does indeed scale like the inverse of the square- 
root of the Schwarzschild radius or like the square-root of the black hole temperature, 
and that the average particle energy (proportional to 7/7/) scales like the square-root 
of the black hole temperature. One important observational effect is that the average 
energy of the outgoing particles is reduced but their number is increased compared 
to direct Hawking emission into vacuum [9, 10]. 



2.6 Photon Emission 

Photons observed far away from the black hole primarily come from one of two 
sources. Either they are emitted directly in the form of a boosted black-body spec- 
trum, or they arise from neutral pion decay. We will consider each of these in turn. 

2.6.1 Direct photons 

Photons emitted directly have a Planck distribution in the local rest frame of the 
fluid. The phase space density is 

The energy appearing here is related to the energy as measured in the rest frame of 
the black hole and to the angle of emission relative to the radial vector by 

E' = 7/(1 — Vf cos9)E . (2.29) 
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No photons will emerge if the angle is greater than n/2. Therefore the instantaneous 
distribution is 

^ _^£>^ ln / 1 _ e - E /2 1/ TA (2 30) 

where the second equality holds in the limit 7/ 3> 1. This limit is actually well 
satisfied for us and is used henceforth. 

The instantaneous spectrum can be integrated over the remaining lifetime of the 
black hole straightforwardly. The radius and boost are both known in terms of the 
Hawking temperature Th, and the time evolution of the latter is simply obtained 
from solving Eq. (2.3). For a black hole that disappears at time t = we have 

1 / 2 \ i/3 

T »v-l(Z) ■ (2 ' 31) 

Here a h is approximately 0.0044 for T H > 100 GeV which includes the contribution 
from gravitons and neutrinos. Starting with a black hole whose temperature is T we 
obtain the spectrum 

i_K = 3 m4 '" '311 £ L p****" ix x * e - . (2 . 32 ) 

dE 7T df V 2u s d f J E ^i n ^ 

Here we have ignored the small numerical difference between a e /^ and a^. In the high 
energy limit, namely, when E ^> 2 / ~ff(To)Tf, the summation yields the pure number 
4(27r 6 /315). Note the power-law behavior E~ 4 . This has important observational 
consequences. 

2.6.2 7T° decay photons 

The neutral pion decays almost entirely into two photons: n° — > 77. In the rest 
frame of the pion the single photon Lorentz invariant distribution is 

d 3 p nm n 
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which is normalized to 2. This must be folded with the distribution of n° to obtain 
the total invariant photon distribution 

E^L = r iE.pi.-L-S ( EE - - PP - - HL) . (2.34) 
d 6 pdt Jm w dE w dt irrriir \ m n 2 / 

After integrating over angles we get 

2 r (2.35) 
p w dE w dt 



d 2 N* „ r °° dE^d 2 N^ 



dEdt 

where E min = (E 2 + m 2 r /A)/E. In the limit E 3> we can approximate E min = E 
and evaluate d 2 N n o/dE 7T dt in the same way as photons. This leads to the relatively 
simple expression 

d _J^_ = ^jh le"^/^/ . (2.36) 
dEdt 7i ^ n 2 K ' 

n=l 

The time-integrated spectrum is computed in the same way as for direct photons 

d£ it df V 2m I^/ / # iin 2 io v ; 

In the high energy limit, namely, when E 3> 2^f(T )Tf, the summation yields the 
pure number 27r 6 /315. 

2.6.3 Instantaneous and integrated photon spectra 

The instantaneous spectra of high energy gamma rays, arising from both di- 
rect emission and from 7r° decay, are plotted in figures 2.10 (for Tf = 140 MeV) 
and 2.11(for Tf = 100 MeV). In each figure there are three curves corresponding to 
Hawking temperatures of 100 GeV, 1 TeV and 10 TeV. The photon spectra are essen- 



tially exponential above a few GeV with inverse slope 2^f(T H )Tf oc ^TfT H . If these 
instantaneous spectra could be measured they would tell us a lot about the equation 
of state, the viscosities, and how energy is processed from first Hawking radiation to 
final observed gamma rays. Even the time evolution of the black hole luminosity and 
temperature could be inferred. 

The time integrated spectra for Tf = 140 MeV are plotted in figure 2.12 for three 
initial temperatures Tq. A black hole with a Hawking temperature of 100 GeV has 
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Figure 2.10: The instantaneous gamma ray spectrum for three different Hawking 
temperatures assuming Tf = 140 MeV. 

5.4 days to live, a black hole with a Hawking temperature of 1 TeV has 7.7 minutes 
to live, and a black hole with a Hawking temperature of 10 TeV has only 1/2 second 
to live. The high energy gamma ray spectra are represented by 

dE ~ 26& ■ (2 - 38) 
It is interesting that the contribution from 7r° decay comprises 20% of the total while 
direct photons contribute the remaining 80%. The E~ A fall-off is the same as that 
obtained by Heckler [9], whereas Halzen et al. [43] and MacGibbon and Carr [8] 
obtained an E~ 3 fall-off on the basis of direct fragmentation of quarks and gluons 
with no fluid flow and no photosphere. 
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Figure 2.11: Same as figure 9 but with Tf = 100 MeV. 

2.7 Observability of Gamma Rays 

The most obvious way to observe the explosion of a microscopic black hole is 
by high energy gamma rays. We consider their contribution to the diffuse gamma 
ray spectrum in Sec. 2.7.1, and in Sec. 2.7.2 we study the systematics of a single 
identifiable explosion. 



2.7.1 Diffuse spectra from the galactic halo 

Suppose that microscopic black holes were distributed about our galaxy in some 
fashion. Unless we were fortunate enough to be close to one so that we could observe 
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Figure 2.12: The time-integrated gamma ray spectrum starting from the indicated 
Hawking temperature. Here Tf = 140 MeV. 

its demise, we would have to rely on their contribution to the diffuse background 
spectrum of high energy gamma rays. 

The flux of photons with energy greater than 1 GeV at Earth can be computed 
from the results of Sec. 2.6 together with the knowledge of the rate density p(x) of 
exploding black holes. It is 

d 3 N Earth _ m\T f f £x p(x) c _ dM/x ^ /(E) (2 3g) 



dEdAdt 26£ 4 J 4vrd 2 (x) 
where d(x) is the distance from the black hole to the Earth. The exponential decay is 
due to absorption of the gamma ray by the black-body radiation [44] . The mean free 
path \ yy (E) is highly energy dependent. It has a minimum of about 1 kpc around 1 
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PeV, and is greater than 10 5 kpc for energies less than 100 TeV. 

We need a model for the rate density of exploding black holes. We shall assume 
they are distributed in the same way as the matter comprising the halo of our galaxy. 
Thus we take 

f>W=fX> x2 + y2 *\ v + Rl (2.40) 

where the galactic plane is the x — y plane, R c is the core radius, and q is a flattening 
parameter. For numerical calculations we shall take the core radius to be 10 kpc. 
The Earth is located a distance Re = 8.5 kpc from the center of the galaxy and lies 
in the galactic plane. Therefore d 2 = (x — Re) 2 + y 2 + z 2 . 

The last remaining quantity is the normalization of the rate density po- This is, 
of course, unknown since no one has ever knowingly observed a black hole explosion. 
The first observational limit was determined by Page and Hawking [45]. They found 
that the local rate density pi OC ai is less than 1 to 10 per cubic parsec per year on the 
basis of diffuse gamma rays with energies on the order of 100 MeV. This limit has 
not been lowered very much during the intervening twenty-five years. For example, 
Wright [46] used EGRET data to search for an anisotropic high-lattitude component 
of diffuse gamma rays in the energy range from 30 MeV to 100 GeV as a signal for 
steady emission of microscopic black holes. He concluded that pi OC ai is less than about 
0.4 per cubic parsec per year. (For an alternative point of view on the data see [47].) 
In our numerical calculations we shall assume a value po = 1 pc~ 3 yr _1 corresponding 
to piocai ~ 0.58 pc~ 3 yr _1 . This makes for easy scaling. Estimating the quantity of 
dark matter in our galaxy as M ha i o /p ha i o = 4.7 x 10 4 kpc 3 means that we could have 
up to 47 x 10 12 microscopic black hole explosions per year in our galaxy. 

Figure 2.13 shows the calculated flux at Earth, multiplied by E 4 . Of course this 
curve would be flat if it were not for absorption on the microwave background radi- 
ation. There is a relative suppression of three orders of magnitude centered between 
10 15 and 10 16 eV. This means that it is unlikely to observe exploding black holes in 
the gamma ray spectrum above 10 14 eV. Even below that energy it is unlikely because 
they have not been observed at energies on the order of 100 MeV, and the spectrum 
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Figure 2.13: The flux of diffuse gamma rays coming from our galactic halo. The 
normalization is po — 1 pc~ 3 yr -1 . The halo is assumed to be spherically symmetric, 
q — 1; the results for a flattened halo with q = 2 are very similar. 

falls faster than the primary cosmic ray spectrum oc E 2J . The curve displayed in 
figure 2.13 assumes a spherical halo, q — 1, but there is hardly any difference when 
the halo is flattened to q — 2. 

2.7.2 Point source systematics 

Given the unfavorable situation for observing the effects of exploding microscopic 
black holes on the diffuse gamma ray spectrum, we now turn to the consequences for 
observing one directly. How far away could one be seen? Let us call that distance d max . 
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We assume that <i max < A 77 for simplicity, although that assumption can be relaxed if 
necessary. Let A^t denote the effective area of the detector that can measure gamma 
rays with energies equal to or greater than E min . The average number of gamma rays 
detected from a single explosion a distance d max away is 

/N(F^F \\ Adet f°° dN ^F- Adct (9A-\\ 
(N,(E > £ min )> - _- I —dE - 4^^- • (2-41) 



max "mm 



Obviously we should have E min as small as possible to get the largest number, but it 
cannot be so small that the simple E~ A behavior of the emission spectrum is invalid. 
See figure 2.12. 

A rough approximation to the number distribution of detected gamma rays is a 
Poisson distribution. 

P(AT 7 ) = i^l e -(^) (2.42) 

The exact form of the number distribution is not so important. What is important is 
that when {N y (E > E min )) > 1 we should expect to see multiple gamma rays coming 
from the same point in the sky. Labeling these gamma rays according to the order in 
which they arrive, 1, 2, 3, etc. we would expect their energies to increase with time: 
Ei < Ei < E 3 < • ■ • . Such an observation would be remarkable, possibly unique, 
because astrophysical sources normally cool at late times. This would directly reflect 
the increasing Hawking temperature as the black hole explodes and disappears. 

It is interesting to know how the average gamma ray energy increases with time. 
Using Eqs. (2.30) and (2.36) we compute the average energy of direct photons to be 
47jTjC(4)/C(3) and the average energy of n° decay photons to be one-half that. The 
ratio of direct to decay photons turns out to be n. Therefore the average gamma 
ray energy is 3.177/(t)7/. This average is plotted in figure 2.14 for 10 5 > t > 10~ 5 
seconds. The average gamma ray energy ranges from about 4 to 160 GeV. 

The maximum distance can now be computed. Using some characteristic numbers 
we find 

. ,„ / -4d« /10GeV\ 3 ' 2 
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Figure 2.14: The average gamma ray energy as a function of the remaining lifetime 
of the black hole. The times spanned correspond to approximately 400 GeV <T H < 
200 TeV. 

If we take the local rate density of explosions to be 0.4 pc~ 3 yr _1 then within 150 
pc of Earth there would be 5 x 10 6 explosions per year. These would be distributed 
isotropically in the sky. Still, it suggests that the direct observation of exploding black 
holes is feasible if they are near to the inferred upper limit to their abundance in our 
neighborhood. We should point out that a search for 1 s bursts of ultrahigh energy 
gamma rays from point sources by CYGNUS has placed an upper limit of 8.5 x 10 5 
p C -3 y r -i |4gj_ However, as we have seen in figure 2.14 and elsewhere, this is not 
what should be expected if our calculations bear any resemblance to reality. Rather 
than a burst, the luminosity and average gamma ray energy increase monotonically 
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over a long period of time. 

2.8 Neutrino Emission 

In this section we focus on high energy neutrino emission from black holes with 
Hawking temperatures greater than 100 GeV and corresponding masses less than 
10 s kg. It is at these and higher temperatures that new physics will arise. Such a 
study is especially important in the context of high energy neutrino detectors under 
construction or planned for the future. Previous notable studies in this area have 
been carried out by MacGibbon and Webber [7, 49] and by Halzen, Keszthelyi and 
Zas [50], who calculated the instantaneous and time-integrated spectra of neutrinos 
arising from the decay of quark and gluon jets. 

The source of neutrinos in the viscous fluid picture is quite varied. Neutrinos 
should stay in thermal equilibrium, along with all other elementary particles, when 
the local temperature is above 100 GeV. The reason is that at energies above the 
electro weak scale of 100 GeV neutrinos should have interaction cross sections similar 
to those of all other particles. Thus the neutrino-sphere, where the neutrinos decouple, 
ought to exist where the local temperature falls below 100 GeV. The spectra of these 
direct neutrinos are calculated in Sec. 2.8.1. 

Neutrinos also come from decays involving pions and muons. The relevant pro- 
cesses are (i) a thermal pion decays into a muon and muon-neutrino, followed by the 
muon decay /i — > ez/ e z/ M , and (ii) a thermal muon decays in the same way. The spectra 
of neutrinos arising from pion decay are calculated in Sec. 2.8.2 while those arising 
from direct or indirect muon decay are calculated in Sec. 2.8.3. 

The spectra from all of these sources are compared graphically in Sec. 2.9. We also 
compare with the spectra of neutrinos emitted directly as Hawking radiation with- 
out any subsequent interactions. The main result is that the time-integrated direct 
Hawking spectrum falls at high energy as E~ 3 whereas the time-integrated neutrino 
spectrum coming from a fluid or from pion and muon decays all fall as E~ A . Thus 
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the fluid picture predicts more neutrinos at lower energies than the direct Hawking 
emission picture. If a microscopic black hole is near enough the instantaneous spec- 
trum could be measured, and its shape and magnitude would provide information on 
the number of degrees of freedom in the nature on mass scales exceeding 100 GeV. 

2.8.1 Directly Emitted Neutrinos 

In this section we first review the emission of neutrinos by the Hawking mechanism 
unmodified by any rescattering. Then we estimate the spectra of neutrinos which 
rescatter in the hot matter. The last scattering surface should be represented approx- 
imately by that radius where the temperature has dropped to 100 GeV. The reason 
is that neutrinos with energies much higher than that have elastic and inelastic cross 
sections that are comparable to the cross sections of quarks, gluons, electrons, muons, 
and tau leptons. Much below that energy the relevant cross sections are greatly sup- 
pressed by the mass of the exchanged vector bosons, the W and Z. Furthermore, the 
electroweak symmetry is broken below temperatures of this order, making it natural 
to place the last scattering surface there. A much better treatment would require 
the solution of transport equations for the neutrinos, an effort that is perhaps not 
yet justified. All the formulas in this section refer to one flavor of neutrino. Our 
current understanding is that there are three flavors, each available as a particle or 
antiparticle. The sum total of all neutrinos would then be a factor of 6 larger than 
the formulas presented here. 

Direct neutrinos 

The emission of neutrinos by the Hawking mechanism is usually calculated on the 
basis of detailed balance. It involves a thermal flux of neutrinos incident on a black 
hole. The Dirac equation is solved and the absorption coefficient is computed. This 
involves numerical calculations [28, 29, 51]. The number emitted per unit time per 
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unit energy is given by 

d 2 7V, dir r„ 1 



(2.44) 



dEdt 2vr exp(E/T H ) + 1 ' 

where T u is an energy-dependent absorption coefficient. There is no simple analytic 
formula for it. For our purposes it is sufficient to parametrize the numerical results. 
A fair representation is given by 

r -4|(" + e# "r j) + i )- < 2 ' 45 > 

and a fa 1.607. The exact expression has very small amplitude oscillations arising 
from the essentially black disk character of the black hole. The parametrized form 
does not have these oscillations, but otherwise is accurate to within about 5%. 

If there are no new degrees of freedom present in nature, other than those already 
known, then the time dependence of the black hole mass and temperature are easily 
found. The relationship between the time and the temperature [Eq. (2.31)] allows us 
to compute the time- integrated spectrum, starting from the moment when T H = T . 
There is a one-dimensional integral to be done numerically. 

j jydir 27ml r E / T o r 4 / 925 \ 

' dx 0.075 + — ^ (2.46) 



dE lQ(47r) 6 a h E 3 Jo exp(x) + 1 \ exp(5 - ax) + l j 

In the high energy limit, meaning E ^> To, the upper limit can be taken to infinity 
with the result 

dN? 29.9ml ^ ml 
dE {4n) e a h E 3 ~ 565£ 3 ' 1 ' ' 

Direct neutrinos from an expanding fluid 

Neutrinos emitted from the decoupling surface have a Fermi distribution in the local 
rest frame of the fluid. The phase space density is 

= <2 - 48) 
The decoupling temperature of neutrinos is denoted by T v . The energy appearing 
here is related to the energy as measured in the rest frame of the black hole and to 



the angle of emission relative to the radial vector by 
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E' = 7„(1 — v v cos 9)E . (2.49) 

No neutrinos will emerge if the angle is greater than ir/2. Therefore the instantaneous 
distribution is 

J2 7yfluid / p2 \ r \ 2rp oo (_i\n+l 

^kr-^Avi) / d(cose)co S ef(E,cose) = r -^Ej: [ - JL - 

dEdt \2it 2 J Jo TiUy ^ n 

(fl - exp[-n£( 7 , - Ul/ )/r v ] + -^L-expl-nE^/Tv]) , (2.50) 

where r v is the radius of the decoupling suface and u = try. Integration over the 
energy gives the luminosity (per neutrino). 

We need to know how the radius and radial flow velocity at neutrino decoupling 
depend on the Hawking temperature or, equivalently, the black hole mass; we have 
already argued that the neutrino temperature at decoupling is about T u = 100 GeV. 
We also know the relation between r v and u v from Eq. (2.20): 

^~f^ 3 . (2.5D 



AttTh \usJ 

The final piece of information is to recognize that each type of neutrino will contribute 
7/8 effective bosonic degree of freedom to the total number of effective bosonic degrees 
of freedom in the energy density of the fluid. After integrating the instantaneous 
neutrino energy spectrum to obtain its luminosity we equate it with the appropriate 
fraction of the total luminosity. 

t fluid _ 7 7r 3 r^X^ 3 — V v _ 7/8 Hawking (2 col 

8 90 (l-^) 3 7 4 106.75 1 ; 

where 

L Hawking = 647r 2 a e// T 2 _ ^53) 

The number 106.75 counts all effective bosonic degrees of freedom in the standard 
model excepting gravitons. Here a e /^ does not include the contribution from gravi- 
tons. This results in an equation which determines v, equivalently u, in terms of the 
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Hawking temperature: 

(3-^W _ 2 i3^V// m6 (Th\ A _ c (Th\ A (2 54) 

(l-^) 3 " 2 427** MrJ - C {tJ ■ (2 ' 54) 

Numerically the constant C — 1.14 x 10~ 5 . 

We are interested in black holes with T# > 71, = 100 GeV. The corresponding 
range of neutrino-sphere flow velocities corresponds to 0.2 < u v < oo. There is no 
simple analytical expression for the solution to Eq. (2.54) for this wide range of the 
variable. At asymptotically high temperatures the left side has the limit 16m®. For 
intermediate values the left side may be approximated by 22w£. We approximate the 
left side by the former for u v > 11/8 and by the latter for 0.2 < u v < 11/8. Thus 

Uv = {22^) fOT Uv < 11/8 (2 ' 55) 
u» = (j^lj for ^ > 11/8. (2.56) 

This approximation is valid to better than 20% within the range mentioned. 

The time-integrated spectrum can be calculated on the basis of Eqs. (2.50), (2.31), 
(2.51) and either (2.54), which is exact, or with the approximation of (2.55-56) which 
results in 



£(-!) 



n+l 



dE 2Va h u 6 Sn , 

C \ 5/7 [T* dT H exp[-nf (71 - Ul )] 



x 



+ 



+ 



\ 5/T [T* dT H exp[-n^( 7 i 



\\22T}) Jt t%> 

T u ( C \ 4/7 [T* dT H /exp[-n^7i] exp[-n^(7i - m)} 



E \22Tt) h 



E \22T*J Jt Tp 7 \ n 2 n 2 



f C V /8 r 
\ 16T, 4 J h 



dT H exp[-n^(7 2 - u 2 )\ 



J Jt„ tV 2 n 

ti 

4/8 



| T v ( C V r°° dT H / exp[-nA 72 ] exp[-n|( 72 -ii 2 )] ] 
E \16T* J Jt+ In 2 n 2 J y 



where T* = 2(^) 2 T U / C 1 ^ . What is interesting here is the high energy limit, E 3> 
10T^ = 1 TeV. This is determined wholly by the first of the two exponentials in Eq. 
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(2.50) with the coefficient of 1. Using also Eqs. (2.31), (2.51) and (2.55) we find the 
limit 

c^ fluid ^ 4185C(6)C 1 / 4 m\T v ^ m 2 T u 
dE ~* 854tt 5 E 4 ~ 1040E 4 ' [ ' ' 

Here we have ignored the small numerical differences between a e J f and ah- This E~ 4 
spectrum is characteristic of all neutrino sources in the viscous fluid description of 
the microscopic black hole wind. 

2.8.2 Neutrinos from Pion Decay 

Muon-type neutrinos will come from the decays of charged pions, namely 7r + — > 
/i + z/ M and 7r~ — > yTv^. In what follows we calculate the spectrum of muon neutrinos; 
the spectrum for muon anti-neutrinos is of course identical. 

In the rest frame of the decaying pion the muon and the neutrino have momentum 
q determined by energy and momentum conservation. 

m 7T = q + yjml + q 2 (2.59) 

Numerically q = 0.2134m,,-. In the pion rest frame the spectrum for the neutrino, 
normalized to one, is 

E^^^LZA. ( 2 . 60 ) 
d 6 p Anq 

The Lorentz invariant rate of emission is obtained by folding together the spectrum 
with the rate of emission of 7r + . 



E- 



<a = r dE „ j-f^-p-p--,) (2.61) 

at JniK \ dE n dt I Anq \ m v J 



d 3 pdt 

The spectrum of pions is computed in the same way as the spectrum of direct neutrinos 
from the expanding fluid in Sec. 2.8.1 but with one difference and one simplification. 
The difference is that the pion has a Bose distribution as opposed to the Fermi 
distribution of the neutrino. The simplification is that the fluid at pion decoupling 
has a highly relativistic flow velocity with Uf « 7/ >> 1. The analog to Eq. (2.50) is 

d 2 N n r 2 / 



dE n dt -rcyj 



&L\n{l-e- E ^f T f) . (2.62) 
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The instantaneous energy spectrum of the neutrino is reduced to a single integral. 

d 2 W v m n T f rj 



^ dE n ln(l-e- E ^f T f) (2.63) 



dEdt 2irq r Yf 

Here E min = m^iE 2 + q 2 )/2Eq is the minimum pion energy that will produce a 
neutrino with energy E. We are interested only in high energy neutrinos with E 
m n , in which case -E^m = m n E/2q is an excellent approximation. The integral can 
also be expressed as an infinite summation 

d 2 N; m.rjTf ~ 1 / nm„E \ 

S^T^J (2 ' M) 

When E ^> 7/7/ only the first term in the summation in Eq. (2.64) is important. 

The instantaneous spectrum is easily integrated over time because the flow velocity 
at pion decoupling is highly relativistic. The spectrum arising from the last moments 
when the Hawking temperature exceeds T is 

dNl 64«fg 3 / 45vrc4 // \ 5/4 m\T s ~ 1 f m^E/4 lf (T )T fq 

' 1 ! zl~I d.r.r 'v " ' . (2.65) 



dE n % a h ml \ 2v? s d f J E 4 ^ n 2 Jo 

In the limit that E 3> 7/(T )T/ the upper limit on the integral may be taken to 
infinity. After ignoring the small numerical difference between a e h ^ and a^, we can 
express the result in terms of the constant C as: 

dNl ^ 464vr / j_X 3 /203C;\ l " m|T> ^ mfj) 
dE 427 VmJ ^ 2d/ j E A 300£ 4 ' 1 ' ' 

This is much smaller than the direct neutrino emission from the fluid because the 
neutrino decoupling temperature T v is much greater than the pion decoupling tem- 
perature Tf. Here a e ^ does not include gravitons and neutrinos. 

2.8.3 Neutrinos from Muon Decay 

Muons can be emitted directly or indirectly by the weak decay of pions: n~ — > /i~z/ e z/ M 
plus the charge conjugated decay. Both sources contribute to the neutrino spectrum. 
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The invariant distribution of the electron neutrino (to be specific) in the rest frame 
of the muon is 



E 



d 3 p' 



(3m M - AE')E' 



(2.67) 



where the electron mass has been neglected in comparison to the muon mass. This 
distribution is used in place of the delta-function distribution of Eq. (2.60) for both 
direct and indirect muons, as evaluated in the following two subsections. 



Neutrinos from direct muons 

The instantaneous spectrum of v e , P e , or arising from muons in thermal equi- 
librium until the decoupling temperature of 7) can be computed by folding together 
the spectrum of muons together with the decay spectrum of neutrinos in the same 
way as Eq. (2.61) was obtained. Using Eqs. (2.62) and (2.67) results in 
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(2.68) 



nE 7/ 7)_ 

where Ei is the exponential-integral function. In the high energy limit, defined here 
by E ^> 27/T/ , the spectrum simplifies to 



d 2 Nf r ^ _ 2r)T f E 



E 
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(2.69) 



dEdt 37T7/ y 27/7) J 
The time-integrated spectrum can be calculated in a fashion analogous to that 
followed in Sec. 2.8.2. Thus 
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In the high energy limit this simplifies to 
dN? T » 11745 
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Neutrinos from indirect muons 



The spectrum of neutrinos coming from the decay of muons which themselves came 
from the decay of pions proceeds exactly as in the previous subsection, but with 
the replacement of the direct muon energy spectrum by the indirect muon energy 
spectrum. In the rest frame of the pion the muon distribution is 



E - 11 



2 — 6 (P - Q) , 



(2.72) 



d 3 p' Anq 2 

which is the finite mass version of Eq. (2.60). Folding together this spectrum with 
the spectrum of pions (2.62) yields the spectrum of indirect muons. 
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This spectrum is now folded with the decay distribution of neutrinos to obtain 
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where 
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The high energy limit E 3> 7/7/ is 
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The integral over time can be done in the usual way. The resulting expression is 
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where Xo± = £±(^o)- The high energy limit is simple and of the familiar form l/E 4 . 

(2.78) 
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2.9 Comparison of Neutrino Sources 

In this section we compare the different sources of neutrinos that were computed in 
the previous sections. All the figures presented display one type of neutrino or anti- 
neutrino. That type should be clear from the context. For example, equal numbers 
of u e , P e , u T , and v T are produced as Hawking radiation and by direct emission 

by the fluid at the neutrino decoupling temperature T v . Only electron and muon type 
neutrinos are produced by muon decay, and only muon type neutrinos by pion decay. 
These differences could help to distinguish the decay of a microscopic black hole from 
other neutrino sources if one happened to be within a detectable distance. 

The instantaneous spectra are displayed in figure 2.15 for a Hawking temperature 
of 1 TeV corresponding to a black hole mass of 10 kg and a lifetime of 7.7 minutes. 
The instantaneous spectra for a Hawking temperature of 10 TeV corresponding to 
a black hole mass of 10 6 kg and a lifetime of 0.5 seconds are displayed in figure 
2.16. There are several important features of these spectra. One feature is that the 
spectrum of direct neutrinos emitted by the fluid, at the decoupling temperature of 
T v = 100 GeV, peaks at a lower energy than the spectrum of neutrinos that would 
be emitted directly as Hawking radiation. The peaks are located approximately at 
\fT v T H for the fluid and at 5T H for the Hawking neutrinos. The reason is that the 
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Figure 2.15: The instantaneous direct neutrino spectra emerging from the fluid with 
neutrino-sphere located at T u = 100 GeV compared to the direct Hawking radiation. 
Also shown are neutrinos arising from pion, direct muon and indirect muon decays 
at a decoupling temperature of Tf = 100 MeV. Here the black hole temperature is 
Tjj = 1 TeV. All curves are for one flavor of neutrino or anti-neutrino. 



viscous flow degrades the average energy of particles composing the fluid, but the 
number of particles is greater as a consequence of energy conservation. In the viscous 
fluid picture of the black hole explosion direct neutrinos are emitted as Hawking 
radiation without any rescattering when Tjj < T u , whereas when Tu > T v they are 
assumed to rescatter and then be emitted from the neutrino-sphere located at T v . 
It is incorrect to add the two curves shown in these figures. In reality, of course, it 
would be better to use neutrino transport equations to describe what happens when 
Tu ps T v . That is beyond the scope of this thesis, and perhaps worth doing only if 
and when there is observational evidence for microscopic black holes. 
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Figure 2.16: Same as figure 2.15 but with Tjj = 1 TeV. 
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Here we would like to investigate the discontinuity in the transition from Hawking 
radiation to neutrino emission from fluid. In order to see the discontinuity, we plot 
the instananeous direct neutrino spectra emitted by the fluid and Hawking radiation 
in figures 2.17, 2.18 and 2.19 for black hole temperatures of 100, 110 and 120 GeV. 
As one can see, the difference in the slope of both process at high energy reaches a 
minimum when the black hole temperature is around Tjj = 110 GeV (figure 2.18). 
We also calculate the total number of neutrinos per unit time for these three graphs 
numerically. For Th = 100 GeV 
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Figure 2.17: The instantaneous direct neutrino spectra emerging from the fluid with 
neutrino-sphere located at T u = 100 GeV compared to the direct Hawking radiation. 
Both curves are for one flavor of neutrino or anti-neutrino. 



For T H = 110 GeV 



For T H = 120 GeV 
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In this energy range the total number of emitted neutrinos from both process is 
roughly the same. Therefore if the transition from Hawking radiation to fluid emission 
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Figure 2.18: Same as figure 2.17 but with T H = 110 GeV. 

happens when T# is around 100 GeV the discontinuity will be small, otherwise the 
discontinuity is big. 

Another feature of figures 2.15 and 2.16 to note is that the average energy of 
neutrinos arising from pion and muon decay is much less than that of directly emitted 
neutrinos. Again, the culprit is viscous fluid flow degrading the average energy of 
particles, here the pion and muon, until the time of their decoupling at Tf ~ 100 
MeV. On the other hand their number is greatly increased on account of energy 
conservation. The average energies of the neutrinos are somewhat less than their 
parent pions and muons because energy must be shared among the decay products. 
The spectrum of muon-neutrinos coming from the decay tt — > \w is the softest because 
the pion and muon masses are very close, leaving very little energy for the neutrino. 




The time-integrated spectra, starting at the moment when the Hawking temper- 
ature is 1 and 10 TeV, are shown in figures 2.20 and 2.21, respectively. The relative 
magnitudes and average energies reflect the trends seen in figures 2.15 and 2.16. At 
high energy the Hawking spectrum is proportional to E~ 3 while all the others are 
proportional to E~ 4 , as was already pointed out in the previous sections. Obviously 
the greatest number of neutrinos by far are emitted at energies less than 100 GeV. 
The basic reason is that only about 5% of the total luminosity of the black hole is 
emitted directly as neutrinos. About 32% goes into neutrinos coming from pion and 
muon decay, about 24% goes into photons, with most of the remainder going into 
electrons and positrons. 
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Figure 2.20: The time integrated neutrino spectra emerging from a microscopic black 
hole. Here the calculation begins when the black hole temperature is T = 1 TeV. 
Either the direct Hawking radiation of neutrinos or the direct neutrino emission from 
a neutrino-sphere at a temperature of 100 GeV should be used. All curves are for one 
species of neutrino or anti-neutrino. 

2.10 Observability of the Neutrino Flux 

We now turn to the possibility of observing neutrinos from a microscopic black hole 
directly. Obviously this depends on a number of factors, such as the distance to the 
black hole, the size of the neutrino detector, the efficiency of detecting neutrinos as a 
function of neutrino type and energy, how long the detector looks at the black hole 
before it is gone, and so on. 

For the sake of discussion, let us assume that one is interested in neutrinos with 
energy greater than 10 GeV and that the observational time is the last 7.7 minutes 
of the black hole's existence when its Hawking temperature is 1 TeV and above. As 
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Figure 2.21: Same as figure 2.20 but with T = 10 TeV. 

may be seen from figure 2.20, most of the neutrinos will come from the decay of 
directly emitted muons. Integration of Eq. (2.71) from E min = 10 GeV to infinity, 
and multiplying by 4 to account for both electron and muon type neutrinos and 
anti-neutrinos, results in the total number of 

4 m 2 T> 



x 10 31 . 



(2.82) 



750 E 3 min 

This does not take into account neutrinos directly emitted from the fluid. For E m[n = 
1 TeV, for example, Eq. (2.58) should be used in place of Eq. (2.71) (see figure 2.20), 
and taking into account tau-type neutrinos too then yields a total number of about 
3 x 10 28 . For an exploding black hole located a distance d from Earth the number of 
neutrinos per unit area is 

2 



N U (E > 10 GeV) = 6700 (^p) km" 2 



(2.83) 
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N V {E > 1 TeV) = 2.5 (^p) km' 2 . (2.84) 

Although the latter luminosity is smaller by three orders of magnitude, it has two 
advantages. First, 1/3 of that luminosity comes from tau-type neutrinos. Unlike 
electron and muon-type neutrinos, the tau-type is not produced by the decays of 
pions produced by interactions of high energy cosmic rays with matter or with the 
microwave background radiation. Hence it would seem to be a much more charac- 
teristic signal of exploding black holes than any other cosmic source (assuming no 
oscillations between the tau-type and the other two species). Second, the tau-type 
neutrinos come from near the neutrino-sphere, thus probing physics at a temperature 
of order 100 GeV much more directly than the other types of neutrinos. 

What is the local rate density pi oca i of exploding black holes? This is, of course, 
unknown since no one has ever knowingly observed a black hole explosion. The first 
observational limit was determined by Page and Hawking [45]. They found that the 
local rate density is less than 1 to 10 per cubic parsec per year on the basis of diffuse 
gamma rays with energies on the order of 100 MeV. This limit has not been lowered 
very much during the intervening twenty-five years. For example, Wright [46] used 
EGRET data to search for an anisotropic high-lattitude component of diffuse gamma 
rays in the energy range from 30 MeV to 100 GeV as a signal for steady emission 
of microscopic black holes. He concluded that pi oca i is less than about 0.4 per cubic 
parsec per year. If the actual rate density is anything close to these upper limits the 
frequency of a high energy neutrino detector seeing a black hole explosion ought to 
be around one per year. 



Chapter 3 



Cosmic Shells 



In this chapter we investigate numerical solutions to the combined field equations 
of gravity and a scalar field with a potential which has two non-degenerate minima. 
The absolute minimum of the potential is the true vacuum and the other minimum 
is the false vacuum. This potential is shown in figure 3.1. The true vacuum of the 
potential is at / 2 and the false vacuum is at f\. The two minima are separated by a 
barrier. 

3.1 Shell Solutions in Static Coordinates 

Consider a scalar field coupled to gravity with the Lagrangian 

c = Y^G R+ l M '"- V[(t)] (3 ' 1} 

where the potential V[4>] has two minima at f\ and ji separated by a barrier. See 
figure 3.1. 

We look for spherically symmetric configurations in which the metric of space-time 
is written as 

ds 2 = --dt 2 + + r 2 (d6 2 + sin 2 9dip 2 ) . (3.2) 
p H 

The functions 0, H, and p depend only on r and t. A tetrad basis is chosen, in a 
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Figure 3.1: The potential V[<j>\. 



region H > 0, as 



e = at , ei = — = ar , e 2 = rao 1 , e 3 = r sin 



(3.3) 

p VH ~ " 

The components of the energy-momentum tensor in the tetrad basis, T ab = e^e^T^, 
are 

Too = ^0 2 + W 2 )+V-[0] 



T 2 2 — T S3 



2\H 



01 



Here dots and primes indicate t and r derivatives, respectively. 
The scalar field satisfies 

'r 2 H 



P- 



d / p ■ \ p d /i 



dt\H J r 2 dr\ p 



' )+V'[<l>]=0 



(3.4) 



(3.5) 
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We introduce the integrated mass function M(t, r) by 



H=l-™i. (3.6) 
r 

The Einstein equations are 

M = I" 4irr 2 drT 00 , (3.7) 
Jo 

^ = _ 4lrGr {|L^ + ^}, (3. 8 ) 

— = -87rGr00' , (3.9) 

§{£(ft+l(?-^)}^-™(^-**-)- <-) 

One of the equations, Eq. (3.10), is redundant as it follows from Eqs. (3.5), (3.7), 
(3.8), and (3.9). 

We shall seek static solutions for which the set of equations reduces to 

0» + r cfr (r)0'(r) = ^'[0] , (3.11) 

M(r) = J\nr 2 dr \^H<f)' 2 + V[<j>]^ , (3.12) 



where 



2 p' H' 

r e fr = \- -jj , 

r p H 

V - = -AitGrcf)' 2 . (3.13) 
p 



Equation (3.11) can be interpreted as an equation for a particle with a coordinate 
and time r. Except for a factor 1/H this particle moves in a potential U[(f>] = —V[<f>]. 
The coefficient r e g(r) represents time (r) dependent friction. 

The potential V[4>] is supposed to have two minima at j\ and We are looking 
for a solution which starts at ~ /i, moves close to / 2 , and comes back to /x at 
r = oo. The particle's potential U[<f>] has two maxima. The particle begins near the 
top of one hill, rolls down into the valley and up the other hill, turns around and rolls 
down and then back up to the top of the original hill. This is impossible in flat space, 
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as r cff is positive-definite so that the particle's energy dissipates and it cannot climb 
back to its starting point. 

In the presence of gravity the situation changes. The non-vanishing energy density 
can make H a decreasing function of r so that r eff becomes negative. The energy lost 
by the particle during the initial rolling down can be regained on the return path by 
negative friction, or thrust. Indeed, this happens. 

Let us set up the problem more precisely. We take a quartic potential with f\ < 
0</ 2 ; 

V'[<j>) = \<j>{<j>-h){<j>-h) , 

V[cf>] = ^(0 - / 2 ){0 3 - \{f 2 + 4/i)0 2 - \f2U2 - 2/i)(0 + h)} ■ (3-14) 

Here V[f 2 ] = 0, and V[0] is a local maximum for the barrier separating the two 
minima. Define / = (|/i| + /2V2 and Af — f 2 — \fi\. In case Af > 0, — /1 
corresponds to a false vacuum with the energy density e = V[fi] = |A/ 3 A/ > 0, 
whereas <p = f'2 corresponds to a true vacuum with a vanishing energy density. As we 
shall discuss in detail below, the positivity of the energy density e plays an important 
role for the presence of shell structure, but the vanishing V[f 2 ] is not essential as we 
see below. In a more general potential it could be that V[f 2 ] > V[/i]. 

We look for solutions with starting at the origin r = very close to f±. There is 
only one parameter to adjust: 0o = 0(0). The behavior of a solution near the origin 
is given by 

= 00 + 2 r 2 + • • • , 02 = \ V'[(f)o\ , 

o 

p — 1 + p 4 r 4 H , p 4 = -47rG02 , 

M = m 3 r 3 + ---, m 3 = —V[(f>o], 

H = l-2Gm 3 r 2 + ■■■ . (3.15) 



Given O the equations determine the behavior of a configuration uniquely. For most 
values of O the corresponding configurations are unacceptable. As r increases, 0(r) 
either approaches (the local maximum of V[4>\) after oscillation, or comes back to 
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cross fi and continues to decrease. Other than the two trivial solutions, corresponding 
to the false and true vacua, we have found a new type of solution. 

There are four parameters in the model, one of which, the gravitational constant 
G = m P 2 , sets the scale. The other three are A, /i, and f 2 or, equivalently, the 
three dimensionless quantities //mp, A///, and A. We have explored only a limited 
region of the parameter space. The moduli space of solutions depends critically on 
//mp and A///. The A dependence can be absorbed by rescaling. The solution 
<f>(r; A) to Eqs. (3.6), (3.11)-(3.13) for a given A is related to the solution for A = 1 by 
0(r; A) = (f)(y/\r; 1). We shall see that nontrivial solutions appear as //mp becomes 
small. 

If 0o < fi, monotonically decreases as the radius increases. On the way H(r) 
goes to zero at a finite radius r m . <f>'(r) may or may not diverge there. To describe 
most general solutions the metric is written as [52] 

do 2 

ds 2 = -Q( p )dt 2 + -^ + R(p) 2 (d9 2 + sin 2 dip 2 ) , (3.16) 

Q{p) 

which is related to Eq. (3.2) by 

R = r , Q = ^ , £ = ±1 . (3.17) 
p z dr \p\ 

If Q(p) tends to zero linearly in p, the space-time has a nondegenerate Killing horizon, 
while if Q(p) tends smoothly to a nonzero value, there is a bag of gold solution studied 
by Volkov et al. [52] in the context of Einstein- Yang-Mills equations. 
As 

\p(r)\ = exp |-47rG dr r0 /2 | , (3.18) 

p can vanish at r m and Q(p) can approach a nonzero value there if <f>' diverges fast 
enough. In this case dR/dp = 0. In other words, as p increases, R(p) reaches a 
maximum R max at p(r m ) and starts to decrease. Eventually R(p) becomes zero at 
Pmax- This would lead to a bag-of-gold solution. It is not clear whether such solutions 
exist in the scalar model under discussion [53]. 

Suppose instead that fi < O < and O is not too close to fi. In the particle 
analogue, the particle starts to roll down the hill under the action of U[4>] = —V[<f>]. 
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It approaches = 0, and oscillates around it. In the meantime H(r) crosses zero. 
We are not interested in this type of solution either. 

Now suppose that 0o is very close to, but still greater than, fa: 0(0) — fi + 50(0) 
with < 5(f)(0)/ f < 1. A schematic of the resulting solution is displayed in figures 
3.2 and 3.3. We divide space into three regions in the static coordinates: region I 
(0 < r < region II (Ri < r < R 2 ), and region III (R 2 < r). It turns out that 0(r) 
varies little from fi in regions I and III so that the equation of motion for may be 
linearized in those regions. In region II the field deviates strongly and the full set of 
nonlinear equations must be solved numerically. This is the region in which we shall 
find shell structure. H(r) deviates from the de Sitter value significantly. In region III 
the space-time is approximately de Sitter again. H(r) crosses zero at r# where p(r) 
remains positive. Consequently Q(p) in Eq. (3.16) also crosses zero linearly so that 
the space-time has a horizon. 
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Figure 3.2: Schematic behavior of 0(r). 
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Figure 3.3: Schematic behavior of H(r) 
In region I the space-time is approximately de Sitter: 



r„ = e , ff=i-£, « = /g^. p = i- (3.19) 

The equation for 0(r) can be linearized with <f>(r) = fi + 5(f>(r). In terms of z = r 2 /a 2 , 

rf 2 /3 5 \ d 1 



^ (1 -^Hi-?Ji-r v r 0=0 ' (3 ' 20) 

where a; 2 = V"[/i]. This is Gauss' hypergeometric equation. The solution which is 
regular at r = is 

6<j>(r) = 50(0) • F(| + iK, § - in, |; , (3.21) 

where 

K =|^|, ^=^(1-^). (3.22) 



63 

We shall soon see that a solution with shell structure appears for wa » 1 with a 
particular choice of 5(f)(0). The ratio of 5<p'(r) to 8(f>(r) is given by 

50(r)"3a 2 V K + 16^(1 + " a 2 l J ' 1 J 

The deviation from /i at the origin, 50(0), needs to be very small for an acceptable 
solution. The behavior of the hypergeometric function for k >> 1 and < z < 1 is 
given by [54] 

F(a + m, a - m, c; z) ~ --7= k^" c z~% + ^(1 - z)^~^~ a exp |2fi;sin" 1 ^fz\ .(3.24) 

The ratio 5<j>(r) /5(f)(0) grows exponentially as r increases like (Ak)~ 1 z^ 1 ^ 2 (1 — z)~ x / 4 
exp |2k sin -1 v^}- At the end of region I, 5(f)/\f\\ needs to be very small for the 
linearization to be valid. The ratio of F'(z) to F(z), J(z) in Eq. (3.23), is given by 

J(z) = . = for k > 1 , < z < 1 . (3.25) 

In region II, <f>(r) varies substantially and the nonlinearity of the equations plays 
an essential role. In this region the equations must be solved numerically. With fine 
tuning of the value of 5(f>(Ri) nontrivial shell solutions will be found. 

The algorithm is the following. First 8(f>(Ri) is chosen and 5(f>'(Ri) is evaluated 
by Eqs. (3.23) and (3.25). To the order in which we work the metric is H(R 1 ) = 
1 — (Ri/a) 2 and p(R\) = 1. With these boundary conditions Eqs. (3.11) and (3.12) 
are numerically solved. 

The behavior of solutions in region II is displayed in figure 3.4. When the specific 
values of the input parameters are chosen to be A = 0.01, f/mp = 0.002, and A/// = 
0.002, then the output parameters are a/l P = 2.365 x 10 7 , k = 3.344 x 10 3 , R\/lp = 
1.76104695 x 10 7 and 50(i?i)/m P = 10" 5 . (Here Z P is the Planck length.) The field 
(f)(r) approaches f\ for r > R 2 . In the numerical integration 5(p(Ri) is kept fixed 
while Ri is varied. Fine tuning to the ninth digit is necessary. If Ri is taken to 
be slightly bigger then 0(r) starts to deviate from /1 in the negative direction as r 
increases. If R\ is taken to be slightly smaller then 0(r) starts to deviate from fi in 
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Figure 3.4: <j>{r) of a solution with f/m P = 1.40 x 10~ 3 , Af/f = 0.002, A = 0.01. 
and r are in the units of m-p and /p, respectively. The maximum value of <fi is smaller 
than f 2 ~ 0.0014 • m P . 

the positive direction, heading for f 2 as r increases. With just the right value the 
space-time becomes nearly de Sitter outside the shell. The value of 5(f) at the origin 
(r = 0) is found from Eq. (3.21) to be 1.40 x 10~ 2440 , which explains why one cannot 
numerically integrate 0(r) starting from r = 0. 

The behavior of the solution in region III is easily inferred. From the numerical 
integration in region II both H 2 = H(R 2 ) and p 2 = p(R 2 ) are determined. In region 
III the metric can be written in the form 



Here M is the mass ascribable to the shell. Once T 00 vanishes H{r) must take this 
form. The value of M may be determined numerically by fitting Eq. (3.26) just 
outside the shell. The location of the horizon, r H , is determined by H(r H ) = 0. The 



H(r) = 1 - 
p(r) = p 2 . 



2GM 



r 



(3.26) 
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field <f>(r) remains very close to /i in region III. 

At this point due caution is required to continue the solution because the static 
coordinates defined in Eq. (3.2) do not cover the whole of space-time. The global 
structure of the space-time and of the solutions is worked out in Sec. 3.4 where 
different coordinates are employed. 

The solutions we found have shell structure at a radius R whose width w is very 
small. The radius R is smaller than the horizon radius r#, but is of the same order 
as a. In the next section we shall see that R becomes smaller as / becomes smaller, 
but remains of order a even in the / — > limit. This implies that the shell structure 
is cosmic in size. 

An estimate of the order of magnitude of w as well as the conditions necessary for 
the existence of the cosmic shell are obtained by a simple argument. Return to Eq. 
(3.11), 

H(f>" + (— + H') <P' + ATvGrH(f)' 3 = V'[<j>\. (3.27) 



V r 

In the shell region (region II, r ~ R) of typical solutions, H is of order one and drops 
sharply. In Eq. (3.27) the 2H/r term is negligible compared with the H' term. The 
remaining terms, H<p", H'<p', and V'[0] are all of the same order of magnitude so that 
f /w 2 ~ A/ 3 or 

"~7xf- (3 ' 28) 

Hence the thickness of the shell is determined by the parameters in the scalar field 
potential. The radius R is smaller than but of the same order as a. The solutions 
exist only if H', which is negative, dominates over AnGrHcf)' 2 . In other words w^ 1 > 
AnGR(f /w) 2 . Making use of R ~ a, a 2 = 9/(16vrG'A/ 3 A/) and Eq. (3.28), one finds 



\mpj a 3 mp 



>30r — • (3.29) 
/ Vmp/ 

We shall see in the next section that these relations are satisfied in the solutions 
obtained numerically. 
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Figure 3.5: H(r) of a solution with f/m P = 1.40 x 10~ 3 , A/// = 0.002, A = 0.01. 
H{r) decreases in two steps in the shell region. H(r) in the de Sitter space (= 
1 — r 2 /a 2 ) is plotted in a dotted line. 

3.2 Numerical Analysis of the Nonlinear Regime 

In the preceding section we solved the linearized field equation for 0(r) in regions 
I and III where the deviation from fi is small, and we sketched the behavior in the 
nonlinear region II where the shell structure appears. In this section we present 
numerical results for region II. As discussed in the preciding section the boundary 
between regions I and II, located at the matching radius Ri, is rather arbitrary, 
subject only to the condition that the linearization is accurate up to that radius. 
Precise tuning is necessary for the pair Ri and 54>(R\) to obtain a solution to all 
the equations. Technically it is easier to keep 5<p(Ri) fixed and adjust the matching 
radius R\. If R\ is chosen too small comes back toward, but cannot reach, f±. What 
happens is that eventually oscillates around <fi = as r increases. If R\ is chosen 
too large comes back to f\ at finite r to roll over it and continues to decrease. 
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Figure 3.6: p(r) of a solution with f/m P = 1.40 x KT 3 , Af/f = 0.002, A = 0.01. p(r) 
decreases in two steps, from 1 to 0.6256, in the shell region. 

For each sufficiently small value of 5<p chosen as the matching value there exists a 
desired solution with R 1 = R c . The critical value R c can be determined numerically 
to arbitrary accuracy. 

One example of such a solution is displayed in figure 3.4 for the parameters 
f/mp = 1.40 x 10- 3 , Af/f = 0.002, A = 0.01 and 5(j)(R 1 ) = 1.0 x 10~ 5 . The 
matching radius R± is fine-tuned to ten digits: Ri/lp ~ 3.7038855228 x 10 7 . The 
shell is very thin compared to the radius of the shell, lying in the region 3.704 x 10 7 < 
r/h < 3.710 x 10 7 . 

In the shell region both H(r) and p(r) decrease in a two-step fashion. See figures 
3.5 and 3.6. Inside and outside the shell H(r) is given by Eqs. (3.19) and (3.26), 
whereas p(r) assumes the constant values 1 and 0.6256, respectively. 

The change is induced by the non-vanishing energy-momentum tensor components 
T 00 = —T22 = —T33 and T n . They are displayed in figures 3.7 and 3.8. The energy 
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Figure 3.7: The energy density Too for a shell solution with //top = 1.40 x 10 -3 , 
Af I f = 0.002, A = 0.01. Too and r are in the units of mp and /p, respectively. The 
energy density is localized in the two shells over the de Sitter background. 

density T 00 has two sharp peaks associated with the rapid variation of 0. The radial 
pressure T n , on the other hand, steps up quickly, remains constant within the shell, 
and then steps down again. The value of Tn is very small (~ 5 x 10 mp) compared 
with the maximum value of Too(~ 2 x 10~ 14 to p ). The contributions of the kinetic 
energy, |/T0' 2 , and potential energy, V [(/>], almost cancel each other. 

In this example 0(r) makes a bounce transition /i — > f 2 — > }\ once. As the 
numerical value of / becomes smaller than 0.002top a new type of solution with 
double bounces /i — > / 2 — > fi — > f 2 — > fi emerge. One example is displayed in figure 



The spectrum of the shell solutions depends on the parameter /. In figure 3.10 
R/a, r H /a, and R/r H are plotted as functions of /. (Finding the solution numerically 
becomes extremely difficult as / becomes smaller.) There are several features to be 



3.9. 
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Figure 3.8: The radial pressure Tu for a shell solution with f/mp = 1.40 x 10 -3 , 
Af I f = 0.002, A = 0.01. Tu and r are in the units of mp and /p, respectively. The 
pressure is higher and constant between the two shells. 

noted. First, with fixed values of A and A/// solutions exist only for / < / max . For 
example, with A = 0.01 and A/// = 0.002 we find that / max ~ 0.006. Second, 
as / goes to 0, r^/a — > 1 and R/a — > 0.8. The relative size R/a of the shell 
remains of order one even in the / — > limit. Third, as / decreases shell solutions 
with multiple bounces become possible. We cannot determine how many times 
can bounce inside the horizon because the numerical evaluation becomes extremely 
difficult when f/m P < 0.001. 

In the potential we are analyzing, <f) = f\ and 4> = fi correspond to the false and 
true vacua, respectively. For the existence of shell solutions the fact V[fi] > V[f 2 ] is 
not crucial, however. In the example displayed in figure 3.4, for instance, <f>(r) swings 
from 0(0) ~ f\ to max < / 2 , and returns to f±. Numerically \^[0 ma x] > V[0(O)]. 
As we are solving classical differential equations, only the form of the potential V[4>] 
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Figure 3.9: A solution with double bounces with f/mp = 1.40 x 10~ 3 ,A/// = 
0.002, A = 0.01. A solution with multiple bounces appears at a smaller shell radius 
than a solution with a single bounce. 



between f\ and max is relevant. The form of the potential for > max does not 
matter. The second minimum can be higher than the first one; V[f\] < V\fy < 
V[4>max\- Before concluding this section we would like to add that we have found no 
solution in which makes a transition from fi to f 2 as r varies from to oo. 



3.3 Global Structure of the Space-Time and Solu- 
tions 

In the preceding section we found novel solutions to a theory with a scalar field 
coupled to gravity in static coordinates. The static coordinates in Eq. (3.2), however, 
do not cover all of space-time. In this section we construct coordinates which allow 
for an extension of the solution to the full space-time manifold. These coordinates are 
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Figure 3.10: The /-dependence of the shell solutions. A = 0.01 and A/// = 0.002 
are fixed. R, r#, and a are the radius of the shell, the horizon length where H(r) 
vanishes, and the horizon length in the de Sitter space, respectively. The ratios of 
various pairs are plotted. Circles, squares, and diamonds correspond to data points 
obtained. 

smooth across the horizon. First we illustrate the construction with de Sitter space, 
and then consider the more general case which is applied to the shell solution. 
The R 1 x S 3 metric of the de Sitter space is given by 



One may regard the de Sitter space as a hypersurface in the five-dimensional Minkowski 
spacetime constrained by the condition 



The metric (3.30) and hypersurface (3.31) cover the entire de Sitter space, whereas 
the static metric covers only half of the space. Furthermore, the static metric has a 



ds 2 = a 2 - 




(3.30) 



yl + vl + vl + vl = a 2 + yl . 



(3.31) 
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coordinate singularity at r = a. 

The relationship among these coordinate systems are easily found. The S 3 metric 
(3.30) can be transformed to a metric conformal to the static Einstein universe, 

„2 



ds 2 
V = 



COS 2 7] 



-drf + dx + sin 2 x 



7T i . 7T 7T 

--2tan- 1 (e^) , - ^ < V < ^ 



(3.32) 



Suppressing S 2 , or 9 and variables, one can map the whole de Sitter space to a 
square region in the x-V coordinates. Null geodesies are given by straight lines at 45 
degree angles. 

The static and hypersurface coordinates are related by 

t 



yo 



= < 



V 'a 2 — r 2 sinh - for r < a 

a 

t 



\Jr 2 — a 2 cosh - for r > a 
a 

Hi = r cos 9 

z/2 = r sin 9 cos <fi 

?/3 = r sin 6* sin 

i 



2/4 



\J d 2 — r 2 cosh - for r < a 
a 

t 



^/r 2 — a 2 sinh — for r > a 

a 



Similarly, Eqs. (3.31) and (3.32) are related by 



(3.33) 



Do = a tan rj 

Hi = a sec r] sin x cos # 

y2 = a sec ?7 sin x sin 9 cos 

yz = a sec 77 sin x sin # sin <fi 

Hi = a sec i] cos x • 



(3.34) 
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The region inside the cosmological horizon (r < a, — oo < t < oo) in the static metric 
corresponds to the left quadrant in the conformal metric (0 < \ < I 71 ") \v\ < 7r /2 — x) 
with the relations 



r sin x 

a cos?] 

t 1 cos x + sin r/ 

- = - In ; 

a 2 cos x — sin 77 



(3.35) 



Similarly, the region outside the cosmological horizon (r > a, —00 < t < 00) in the 
static metric corresponds to the upper quadrant in the conformal metric (0 < x < 11 1 
\v\ > X ~ with the relations 



r sin x 

a cos i] 

t 1 sin 77 + cos x 

- — - m 

a 2 sin 7/ — cos x 



(3.36) 



See figure 3.11. 

For the general static metric given by Eq. (3.2), with H = H(r) and p = p(r), 
the construction of the conformal coordinates analogous to Eq. (3.32) proceeds as 
follows. We suppose that H(r) has a single zero at r H whereas p(r) > 0. The new 
radial coordinate is defined by 



r*(r) = < 



/ dr 

Jo H(r') 



, P(r') 



dr' 

H{r') 



for r < r H 



for r > ru- 



(3.37) 



It has a logarithmic singularity at r#, diverging there as r* ~ — |Mn |r — r#| where 
6 = —2p(r H )/H'(r H ) > 0. New coordinates are introduced by 



tan w 



+e 



+t/b-r*/b 



tant> 



Te -t/b e -r*/b 



(3.38) 
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Figure 3.11: Penrose diagram of the de Sitter space. 



The upper sign is for r < r H and the lower sign is for r > r H . The u-v coordinates 
are related to the U-V coordinates in the Kruskal-Szekeres (KS) coordinate system 
[55] for the Schwarschild solution, and to the Gibbons- Hawking (GH) coordinate 
system [56] for de Sitter space. The connection between the KS and GH coordinate 
systems has already been discussed in [25] . Essentially tan u — U + V and tan v = 
V — U. The metric becomes 

4b 2 F(u,v) 



ds 2 = 



cos 2 (w + V 



-dudv + r 2 dQ 2 



H 

F(u,v) = — -(l — tan u tan i>) (l — cot ucotv 



4p 
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H / , 2r,x „ \r<r H 



(l ± cosh —g-j f° r | 



v v h) ~y r > rH - ( 3 - 39 ) 

The static metric covers the region interior to the bounding lines u — 0, u + v — |7r, 
and u — v = \n. The horizon in the static metric, r = m, corresponds to the 
single point u — v — 0. The function F is non-vanishing and finite there as e 2r *l h ~ 
l/|r — r H \. In the w-f coordinates the metric is regular on u = so that the extension 
to the region u < can be made, whereas the static metric covers only half of the 
space. 

Applied to the de Sitter space we find 

' a — r 



r H = b = a, e' r " /b 
and the metric (3.32) is recovered by 



a + r 



1/2 

, F = l, (3.40) 



7T 

rj = u + v , x = f — -U + — . (3-41) 

The shell solution found in the preceeding sections can be extended to the entire 
space-time. In the region III (R 2 < r) defined in the static metric H(r) is given by 
Eq. (3.26), and the location of the horizon is determined by H(r H ) = 0. In region 
III, 5<f>(r) is very small so that its equation of motion can be linearized. We divide 
region III into two; region Ilia (i?2 < r < 2r# — R2) and region Illb (r > 2r# — R2). 
In region Ilia, H(r) can be approximated by 

H{r)~A(l-£) , (3.42) 



r 2 H 



/ 1 



where 

l/3r f-l). (3.43) 



2 V a 2 

The coefficient A has been chosen such that both Eqs. (3.26) and (3.42) have the 
same slope H'{j-h) at the horizon. In the examples described in Sec. 3.3, errors 
caused by Eq. (3.42) are less than 15% in region Ilia. Now we write the linearized 
version of Eq. (3.11) in terms of y — 1 — (r 2 /r 2 H ), 
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Since 5(j)(r) must be regular at r = ru (y — 0), the solution is 

50( r ) = 5(f>{r H ) • F(| + i«, | - i«, 1; y) . (3.45) 

The normalization 5^{th) must be such that 5<f>(r) matches at r = R 2 with the 
value determined by numerical integration in region II. Essentially S(f>(r) decreases 
exponentially when R 2 < r < rn- 

Near r = R 2 the hypergeometric function behaves as 

F(f + i«, | - i«, 1; y) ~ -4= 2T 1/4 (1 - 1/)" 1/2 exp {2/isin- 1 ^} (3.46) 



so that 

8<ft(r) 2r k 



Wr) rl ^y(l-y)' (3 ' 47) 
The value of the right-hand side of Eq. (3.47) at r = R 2 can be compared with 
the value obtained by direct numerical integration in region II. In one example with 
f/mp = 2 x KT 3 , Af/f = 2 x 1(T 3 , A = 0.01, R 2 = 1.765 x 10 7 , a = 2.365 x 10 7 , r H = 
1.985 x 10 7 , k = 3344, and k = 4481. The numerical value for 5<f)'{R 2 )/5<t){R 2 ) is 
-0.000848, whereas the value from Eq. (3.47) is -0.000985. With the uncertainty 
in the value of k caused by the approximation (3.42) taken into account, one may 
conclude that the agreement is rather good. To gauge the difficulty of determining 
the solution numerically for all values of r, we note that 84>(r H ) ~ 1CT 1850 • 5(f)(R 2 ). 

The most important observation is that the solution 0(r) is regular at the horizon, 
r = rut though the static metric (t,r,0,(f>) is not: r = r# in the static metric is a 
coordinate singularity. As seen above, r = r# corresponds to u — v — in the 
conformal metric (3.39) which is non-singular. Near the horizon tan u tan v ~ r — r# 
so that 

= h + H{r H )[l - ^-(k 2 + A) uv + . . . | . (3.48) 

Furthermore, as tanwtanf = e~ 2r * ( - r ^ b , the static solution <f)(u,v) depends only on 
tan u tan v in the entire space-time. This implies that 

<f>(-u, -v) =<i>(u,v) . (3.49) 
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Figure 3.12: Global structure of the shell solution in the conformal coordinates. The 
mirror image appears in the other hemisphere. 

In the rj-x coordinates (S* 3 metric) defined by Eq. (3.41) 

<f>(v, x) = 0(-?7, 7T - X) = (f>(v, * ~ X) ■ (3-50) 

This relation shows that there is a mirror image of the shell structure in the other 
hemisphere in the S 3 metric. See figure 3.12. The shaded regions in the figure 
represent the shells. The size of the shells is invariant, as is obvious in the static 
coordinates. In the S 3 or conformal metric an observer at the center (x = or 7r) 
sends a light signal toward the shell. Another observer at the shell, upon receiving 
the signal from the center, sent a signal back to the observer at the center. The total 
elapsed time, in terms of the proper time of the observer at the center ar where r is 
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related to rj by Eq. (3.32) at x — 0, is a\n(a + R)/(a — R), independent of when the 
signal is emitted initially, as straightforward manipulations show. 

Finally, in the region outside the horizon in the static metric, or in the upper or 
lower quadrant in the Penrose diagram, 5<p{r) decays as r -3 / 2 . The proof goes as 
follows. For r 3> r H , H ~ 1 — (r 2 /a 2 ), so that 5<J)(r) satisfies Eq. (3.20). Hence it is a 



linear combination of U\ — F(| +ik, 
These can be written as 



IK 



' 2' 



z) and u 2 = z- 1 ' 2 F(\ + 



IK 



' 4 



■IK, 2? z). 



Ui 

U 2 



ir(3 



IK) 



. r(i-m)- 2 

|r(f + m)- 2 
r(i + m)- 2 



0Fr( 



From this one can write 



'_2)(3/4)+»k 



7rr(+2i/c) 

U3/4)-i« 



/x + A^y 72 sin (2k In ^ + 5) 



2m) /q i 1 

i + i/c, l + 2i/c;- 



F(| - m, | - m, 1 - 2m; -) . (3.51) 



(3.52) 



when r th- 



3.4 False Vacuum Decay 

If the region inside the shell is in a false vacuum state one should consider its 
quantum decay to the true vacuum state. The lifetime of the false vacuum may be 
determined semiclassically using the methods of Coleman et al. without [57] or with 
[19] gravity taken into account. The rate per unit volume for making a transition 
from the false vacuum to the true vacuum is expressed as 

£ = Ae- B / n [l + 0(h)] (3.53) 

where Planck's constant has been used here to emphasize the semiclassical nature 
of the tunneling rate. For the potential being used in this paper we find the 0(4), 
Euclidean space, bounce action (neglecting gravity) to be 
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where the radius of the critical size bubble which nucleates the transition is 

For any sensible estimate of the coefficient A the lifetime of the false vacuum will 
exceed the present age of the universe when the condition 

A^ 3 <1 (3.56) 

is satisfied. This calculation is based on the thin wall approximation, which is valid 
when the critical radius is large compared to the coherence length of the potential, 



namely, l/y\V"\. This condition translates into 

A/«6/. (3.57) 

With gravity included the bounce action is 

B = — . . (3.58) 

[1 + (p e /2R)*\ 2 

Gravitational effects are negligible when p c < R. When this condition is not fulfilled 
the shell radius is too small to accommodate even a single nucleation bubble and 
therefore nucleation is further suppressed. 



Chapter 4 



Summary 



The increasing energy of the gamma rays and neutrinos emitted by exploding micro- 
scopic black holes and the disappearance of such point sources in a certain period 
of time are unique characteristics because astrophysical sources normally cool at late 
times. This would directly reflect the increasing Hawking temperature as the black 
hole explodes and disappears. We should emphasize in particular the usefulness of 
distinguishing between the electron and muon-type neutrinos and the tau-type ones. 
The latter are much less likely to be produced by high energy cosmic rays. If the 
rate density of exploding microscopic black holes in our vicinity is anywhere close to 
the current limit based on gamma rays, it should be possible to observe them with 
present and planned large astrophysical neutrino detectors. Observation of high en- 
ergy neutrinos, especially in conjunction with high energy gamma rays, may provide 
a window on physics well beyond the TeV scale. The results of this research will bear 
directly on the observational program on astroparticle physics that involves gamma 
ray detectors such as Milagro and neutrino detectors such as AMANDA and IceCube, 
as well as detectors for ultra high energy cosmic rays, such as Auger, to observe and 
detect possible signals of microscopic black holes. 

Still, there is much work to be done in determining whether the matter surrounding 
a black hole can reach and maintain thermal equilibrium. The previous calculations 
should be improved by including the effect of gravity, which was neglected in the 
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calculation, which may lead to corrections of up to 10%. The equation of state 
should be improved, and the viscosities computed using the relaxation times for self- 
consistency of the transition from viscous fluid flow to free streaming. Also, there 
should be a more fundamental investigation of the relaxation times starting from the 
microscopic interactions. Our next step is to calculate the flux of other particles such 
as antiprotons and positrons radiated by microscopic black holes in order to help the 
observation of these black holes. Positrons and antiprotons are expected to freeze out 
at a temperature around 100 MeV. Another worthwhile project is to carry out cascade 
simulations of the spherically expanding matter around the microscopic black hole at 
a level of sophistication comparable to that of high energy heavy ion collisions. This 
project is much more complicated than the cascade simulation in heavy ion collision, 
though, because we need to deal with a much wider range of energies and particles 
involved in exploding microscopic black holes. 

The study of microscopic black holes might well lead to great advancements in 
fundamental physics. Because the highest temperatures in the universe exist in the 
vicinity of microscopic black holes, matter at extremely high temperatures can be 
studied, and physics beyond the four dimensional Standard Model and above the elec- 
troweak scale can be tested. In addition, because it is believed that baryon number 
is violated at high temperatures, the study of microscopic black holes could possibly 
answer the question of why our universe became matter-dominated. Because micro- 
scopic black holes explode, they are an ideal model for studying the Big Bang and 
the birth of our universe. Finally, the study of microscopic black holes will help us 
to determine whether they are the source of the highest energy cosmic rays. The 
origin of these cosmic rays is still one of the biggest mysteries today. Observation 
and experimental detection of microscopic black holes with satellite and ground based 
neutrino observatories may lead to a novel paradigm for understanding new sources 
of the highest energy cosmic and gamma rays in the Universe. Observation and ex- 
perimental detection of microscopic black holes will be one of the great challenges in 
the new millennium. 
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In the second part of this thesis we have reported the discovery of shell-like so- 
lutions to the combined field equations of gravity and a scalar field with a double- 
minima potential. These solutions exist in a space that is asymptotically de Sitter. 
The range of parameters which allow such solutions are very broad. If anything like 
these structures exist in nature they most likely would have been created in the early 
universe and are therefore cosmological. We know of no other way to produce them. 

To make matters interesting, let us suppose that the cosmological constant sug- 
gested by recent observations of distant Type la supernovae [58] arises from the 
universe actually being in a false vacuum state. A best fit to all cosmological data 
[59] reveals that the present energy density of the universe has the critical value of 
e c = 3Hq/8itG, with one-third of it consisting of ordinary matter and two-thirds of 
it contributed by the cosmological constant. Suppose that the cosmological constant 
arises from a potential of the form we have been analyzing. Then with a present value 
of the Hubble constant of H = 65 km/s-Mpc we find that 

a = = 1.7 x 10 26 m (4.1) 

and so 

(A/ 3 A/) 1/4 = 2.4 x 1(T 3 eV . (4.2) 

This is a constraint on the parameters of the potential, A, / and Af. Although they 
cannot be determined individually from this data, we can place limits on them such 
that a shell structure might arise. Recalling (3.29) we find 

/ < 30A~ 1/6 MeV, 

w > 6.3A~ 1/3 fm. (4.3) 

This is pure speculation of course. A cosmological constant, if it exists, may have 
its origins elsewhere. But if it does arise from a false vacuum, a variety of questions 
immediately present themselves. Is a new field, not present in the standard model 
of particle physics, whose only purpose is this? Where does the energy scale of 2.4 
meV come from? Why should V[<f>] have a global minimum of 0, especially when 
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quantum mechanical fluctuations are taken into account? To these questions we have 
no answers. 



Appendix A 



Numerical Calculations 



In the numerical calculations of Eqs. (2.11), (2.13), (2.14), (2.17) and (2.18), we use 
the unitless variables t = T/Th and x = r/r s . In terms of these unitless variables, we 
can rewrite the equations as 



St 



4vr2 
90 



101.5, t EW < t, 

56.5 + 45 e"^-')/*, t QCD < t < t EW , 



2 + 3.5e' 



-me/T h 



Vt = 



+ 27.25 Q-«QCD-t)/t^ t<t QCD , 
82.5 

St, 



lfolOl.5 



and 



d 



du 
dx 



u 

x 



Tx [x USt] 



J - (4ir)(au | 

8 (du m n 
9^* I dx x , 



du 


2u\ 


dx 


x J 


2 

+ c 


( du 
\ dx 



AitL/Tfj 



2u 

x 



(A.l) 

(A.2) 

(A.3) 
(A.4) 

(A.5) 



where s t = s/T^, rj t = r]/T^ and ^ = £/Tjj. Here L = 647r 2 a fe T| where a h = 
4.43 x 10 -4 includes neutrinos and gravitons. In the case of calculating gamma ray 
spectra, gravitons and neutrinos should be ignored since gravitons never contribute 
to the fluid and neutrinos freeze-out at a temperature around 100 GeV. In the case 
of calculating neutrino spectra, only gravitons need to be ignored for direct neutrinos 
which freeze-out at 100 GeV. For the indirect neutrinos both gravitons and neutrinos 
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need to be ignored. In our calculations these small numerical differences in have 
been neglected. 

In order to find the asymptotic behavior of u and t in the limit where x — > 0, we 
look for a behavior of the form u = Ax a and t = Bx 13 with 7 w 1. For black holes 
with temperatures greater than T EW = 100 GeV, we can write the equation of state 
as 

4 

s t = -at 3 , T) t = b s t 3 and & = b B t 3 , (A. 6) 

where a = ^ 101.5, 65 = 82.5 and b B = b s x 10~ 4 . The best solution in the small x 
limit is 



a 



P = 
A = 



and 



2 

5' 

3 

~5' 

5a J 47rL 

tt(8o 5 + 144o B ) { If 



x -a + 47r(-6 5 — —6 



12, s 5a \ 1 / 5a \ 1 1 1/5 



3~ ' '"'5^ 5 "• B/ tt(86 s + 1446^)/ \n(8b s + 1446 s ) 



f 4vrL /4 „ .4, 12, , 5a 

B = \^r\n a + ^{-bs-—b B )- 



1% V 3 5 5 vr(86 5 + 1446 B ) / 

( 5a \~M 1/5 

X {n(8b s + lUb B ) J } • (A ' 7) 

Typical numerical values for -u and t in this limit are u = 0.105715a; 2 / 5 and t = 
1.31521a;- 3 / 5 . 

In order to find the asymptotic behavior of u and t in the limit where x — > 00, 
we again look for solutions of the form -u = Ax a and t = Bx 13 with 7 ^ m. In this 
limit the equation of state can be written in the same form as in Eqn. (A. 6), where 
a = t^2, 65 = 82.5^^-5 and b B = b s x 10~ 4 . The asymptotic solution in the large x 
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limit is 



1 

Q = 3' 



P 
A 



2 

~3' 

9a J AnL 

7r(326 5 + 4416 B ) \ If 



x -a + 47r(-6 5 - -b B ) 



8, 7, , 9a \-V 9a \ 2 1 1/6 



3 v 9 3 y 7r(326 s + 4416 B ); Vtt(32& s + 441& s ) 



and 



f 4vrL (A / 8 1 7, , 9a 



T| V 3 9 3 7r(32fe 5 + 441fo B ) y 

X U(326 s + 441mJ J ' (A ' 8) 

The typical numerical values for u and t in this limit are u = 0.184460a; 1 / 3 and 
t = 5.09763x- 2 / 3 . 

To solve Eqs. (A.1)-(A.5) we use a 4th order vector Runge-Katta method. In 
our program we apply a shooting type method to find the numerical solution over 
a wide range of x such that u has the correct asymptotic behavior when x — > oo. 
In the program this is done by comparing the numerical solution for u to its correct 
asymptotic behavior at large x. The initial value of t at x — 0.1 is fixed and the initial 
value of u at x = 0.1 is changed by the program until the correct solution is obtained. 
In order to check the validity of the numerical solution, we calculate 4nr 2 T 0r / L using 
finite difference method. The result, which is very close to 1 over a wide range of x, 
is plotted in figure A. 1 for a black hole temperature of 10 TeV. For more information 
see the numerical recipes mbhwind.c and subroutinewind.h attached to this thesis. 
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Figure A.l: 4iir 2 T 0r / L as a function of r/r s for a black hole temperature of 10 TeV. 
The curve begins at r s and terminates when the local temperature reaches T = m e . 
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